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APPLICATIONS 


New Departures in High-Flux 
Research Reactors 


By Bernard |. Spinrad 


The present generation of high-flux research reac- 
tors —heterogeneous, water cooled and moderated 
(light and heavy), with projected beam and irradia- 
tion fluxes of 1 to 5 x 10" thermal neutrons/(cm?) 
(sec) —represents the penultimate in research-re- 
actor performance, according to their designers. 
That is, although they may be improved by design 
refinements as we learn more about actual system 
performance with regard to hot-channel factors, 
coolant flow patterns, burnout safety margins, and 
so on, it is unlikely that this improvement could 
lead to more than a factor of ~ 3 in available flux. (If 
this were not the case, these improvements would 
be the subject of vigorous research and development 
and would be incorporated in the most recent pro- 
posals.) 

To go further in obtaining greater neutron fluxes 
for experimental purposes, it is necessary to over- 
come the limitations of present reactor types and to 
find new reactor concepts that are more amenable 
to very-high-flux design. After discussing methods 
by which present limitations can be overcome by 
relatively moderate design changes, this article con- 
siders a number of specific new concepts. Even 
though several of the new approaches might be able 
to deliver ~10"" neutrons/(cm?)(sec), they, too, have 
foreseeable limits. It is likely that a flux greater 
than 2 x 10" will never be achieved in a steady-state 
research reactor, However, if it were possible to 
“hire” the Jackass Flats hydrogen pumping capability 
for an occasional—and, indeed, very expensive — 
hour, fluxes up to 5 x 10” might be contemplated with 
reactor designs appropriate for 50 Mw/liter opera- 
tion. A similar complex (Fig. 1) that uses stored 
helium to cool a reactor (Fig. 2) might be more at- 
tractive from the standpoint of operating cost. He- 
lium would be safer to handle than hydrogen, and 
there would be a greater variety of reactor ma- 
terials that can tolerate it. 
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Overcoming Limitations 


The major factor that limits the advance of fluxes 
in present types of research reactors is the power 
density, which is restricted by conflicting require- 
menis of heat transfer, heat transport, and core ma- 
terials. Designs now have coolant flows of 6 to 15 
m/sec, core pressure drops of 5 to 10 atm, and core 
power densities of 1 to 3 Mw/liter. Typical heat- 
transfer rates are ~200 watts/cm’, and system 
pressures are tens of atmospheres to suppress boil- 
ing. The materials and heat-transfer technologies do 
not permit much increase in power density with 
present reactor concepts. But, for experiment ir- 
radiation purposes, we can (1) alter the fissile ma- 
terial to obtain greater neutron-production rates, (2) 
take greater advantage of geometry, and (3) use the 
core flux more effectively. These steps lead to con- 
sideration of flux traps, core life, criticality, and 
power shape and their influence on reactor design 
and operation. 


ALTERING THE FISSILE MATERIAL 


Per fission yielding v neutrons, the critical reac- 
tor core absorbs v/k.. neutrons, where k. is the in- 
finite multiplication factor of the core. Therefore the 
excess neutrons produced per fissionare v(k,,—1)/k.. 
Without parasitic absorption, the maximum value k. 
can assume is7. Table 1 lists values of v(7—-1)/n, 
which represents the maximum number of neutrons 
that can be made available per fission, for three 
fissile isotopes. It is apparent that *°*Pu is a very 
superior fuel in the fast-energy range. Uranium-233 
is the best fuel in the intermediate-energy range, 
and *5%py and **%y are about equally good in thermal 
reactors. 

Unfortunately, *“*U is the most available fuel and 
is used preferably. Its use in an intermediate-reac- 
tor core results in only half as many neutrons being 
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Fig. 1 Helium system for cooling reactor (Fig. 2) that might produce 5 x 10'’ neutrons/(cm?) (sec). Coolant 


storage complex might be dubbed “Helium Hill.”’ 
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Fig. 2. Research-reactor concept capable of 50 Mw/liter 
operation to produce 5 x 10!” neutrons/(cm?)(sec). 


available as in a plutonium-fueled fast core, and 
two-thirds as many as in a *y- or **pu-fueled 
thermal reactor. Uranium-233 is the least available 
fuel. It would appear worthwhile, therefore, to ex- 
plore the possibility of using ***Pu, which is plenti- 
ful, to fuel research reactors. The desirability of 
39Pu is enhanced by its high neutron-absorption 
cross section at all energies, which makes it rather 
easy for the infinite multiplication factor to ap- 
proach n. 


Table 1 AVAILABLE NEUTRONS PER FISSION [v(n—1)/y] 








Energy range 233 Ul 239Pu 
Thermal 1.41 1.26 1.43 
Intermediate 1.39 0.96 1.28 
Fast 1.50 1.41 1.90 





TAKING ADVANTAGE OF GEOMETRY 


Experiments are generally located in the reflector 
of the reactor system. To examine the geometrical 
factors, let us take the reactor as a sphere emitting 
Q neutrons, and let 7 be the reflector age, a the core 
radius, L the reflector diffusion length, and D the 
diffusion coefficient for thermal neutrons. To enhance 
the ratio of thermal neutrons to high-energy neu- 
trons, the experiment is usually located in the re- 
flector at a distance of about V7 from the core edge. 
The unperturbed thermal flux is usually within a 
factor of 2 of QL/(a+ V7)(a+L +V7)D. 

The properties of moderators are such that the 
maximum flux is obtained in light water when a =15 
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cm, in beryllium or BeO when 15 cmx<a < 30cm, 
and in D,O when a=30 cm. More importantly, the 
flux increases for each moderator asa decreases— 
i.e., the diffusion processes merely serve to smear 
out the basic fr? law. With moderator optimization 
taken into account, the unperturbed neutron flux 
follows a 1/r law up to about 1 m, whereas current 
(more appropriate for large beam setups) follows 
a 1/r? law. 

Thus it is obvious that the more the source can be 
compressed the more flux we will have—so the 
incentive is for high power density. 

To some extent this argument may be slightly 
softened by the use of “external flux traps” or a 
“neutron hohlraum.” External flux trapping, a means 
of leading more than the geometrical share of neu- 
trons to an experiment, is usually done by putting 
some hydrogenous moderator around the experiment 
(when the reflecting material is beryllium or D,0). 
In a neutron hohlraum the core and a small quantity 
of reflector are placed in a large cavity in a block of 
dense moderator. Even in these cases the advantages 
of high-power-density cores persist. 


USING THE CORE FLUX 

Under circumstances where the backgrounds canbe 
made acceptable, it is desirable to get close to the 
source to obtain high flux for an experiment. In fact, 
in principle, it is possible to extract neutrons di- 
rectly from the core. In this circumstance the pre- 
vious arguments are not valid. Instead, one must 
examine core criticality and power density more 
closely. 

If the core power density is PD (in megawatts per 
liter) and the macroscopic fission cross section is 
=,, then the flux is 


PD 
b= Fe 3.25 x 10% (1) 


The macroscopic fission cross section is related to 
the critical mass and volume by 


=e) () 


where Ny = Avogadro’s number 
o, = the microscopic fission cross section 
A = the molar mass 
M = the critical mass 
V = the core volume 


Then 


g = (3:25 x 1074 55 (5) 


No o f M 


(3.25 x 10") A (P 
eg M (3) 
ove - 
where P is total power. 
Thus, to get high core flux, we attempt to attain a 
highly dilute core (small M/V) or a small critical 
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mass, depending on whether power density or power 
is limiting. Additionally, we may get high epithermal 
fluxes by making the core undermoderated, and 
thereby achieving small values of 0;. However, this 
factor is counterbalanced by the fact that critical 
masses of intermediate reactors tend to be large. 


FLUX TRAPS 

The ultimate in penetrating the core is to use the 
flux at its center. This poses no problem for lightly 
loaded reactors. However, heavily loaded small re- 
actors immersed in a moderator uave peak core flux 
at the outer edge and peak thermal flux in the mod- 
erating reflector at ~1 root age from the core. 

This effect also can be observed in the center of 
the core. Geometrically, if one takes a slab core and 
deforms it into a cylindrical or spherical annulus, 
the thermal-neutron-flux peaks on both sides per- 
sist—and one of these peaks now can be at the cen- 
ter of the annulus. This peak is known as a flux trap 
and has been well investigated.!~® The most definitive 
study is that of Ergen.*° However, it is now well 
known that the performance of a flux trap depends 
strongly on core configuration. 

The following heuristic points have emerged as a 
result of the many investigations of moderator islands 
in research-reactor cores.°® 

e There is a limiting ratio of core-to-moderator 
blackness, below which no flux peak in a moderator 
island is possible. This does not mean that experi- 
mental fluxes in moderator islands embedded in cores 
of low absorption are necessarily low; such cores 
are characterized by already well-thermalized spec- 
tra. The flux trap in such cores has degenerated to 
a simple filter removing virgin-neutron background. 
Small quantities (dimensions small compared to 
thermal-neutron-diffusion length in the moderator) of 
materials such as D,O or beryllium are optimal. In 
fact, maximum experimental flux per unit core power 
density can be obtained from such systems. 

e The flux-trap principle peaks the thermal-neu- 
tron flux most dramatically when cores of high black- 
ness to thermal neutrons are used, which means 
that, in the core, thermal-neutron flux is strongly 
depressed. 


It follows from these two principles that the actual 
value of thermal-neutron flux in the center of a mod- 
erator island only mildly depends on core blackness 
and, in fact, is higher for cores of low blackness: the 
flux produced by thermalization of fast neutrons in 
the island is essentially additive to core flux. 

The preceding argument assumes that “other things 
are equal.” In practice, other things are not equal. 
For example, two idealizations of the moderator 
island assume that (1) the total source strength is 
constant (i.e., fast-neutron current at the edge of the 
island = Q/ 41a’, where Q is the available neutrons) 
and (2) fast current is constant at the edge. It turns 
out that geometrical factors prohibit flux peaking in 
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the former case, which is appropriate to large power- 
limited cores. Per unit source the maximum flux is 
attained with a vanishing island. However, the source 
entering a small island is more nearly described by 
a constant entering current from a power-density - 
limited reactor (Case 2). Nevertheless, this current 
may be made somewhat larger by maximizing the 
inward leakage to the flux trap relative to outward 
leakage. To do this: 

e The core should be as thin as possible. 

e The core should be externally reflected by a good 
fast-neutron reflector, particularly one with a high 
scattering cross section for epithermal neutrons. 
Beryllium and nickel are recommended. 


These two arguments can also be arrived at by 
simply examining the contribution to fast-neutron 
current entering the island from fission sources 
originating at various points in the core. The con- 
tribution is greatest for sources close to the island 
and can be somewhat enhanced by a good external 
reflector. The two points of view are therefore func- 
tionally related. 

So far, we have not commented on steps to maxi- 
mize flux within the island. The following points 
should be considered: 

e The “fast” neutrons leaking in should actually 
be as soft as possible. 

e The neutron-absorption cross section of the flux- 
trap medium should be low, and the diffusion coef- 
ficient for thermal neutrons in it should also be low, 
to permit high flux buildup and to inhibit thermal- 
neutron migration from trap to core. 


All these criteria are difficult to satisfy at once. 
For example, H,O, D,O, and beryllium can all func- 
tion as flux-trap media of approximately equal ef- 
ficiency, but the dimensions of the optimum trap 
vary greatly and the optima emphasize different 
properties of the media (the low age and thermal- 
neutron-diffusion coefficient of water, the low ab- 
sorption by D,O, and the intermediate properties of 
beryllium). 

Although, ideally, the flux trap is completely sur- 
rounded by core, in practice not too much harm is 
done by penetrating the core with a beam tube. As 
with irradiation experiments, not too many neutrons 
must be extracted by means of a beam or the flux 
will be depressed. About half the neutrons in the trap 
may be used. 


CORE LIFE AND REFUELING 


Very-high-flux research reactors will operate at 
quite high power—100 to 1000 Mw, and maybe even 
more. At these powers, fissile fuel is destroyed quite 
rapidly. At 100 Mw, about 1 kg is destroyed in 10 
days, and at 1000 Mv, in 1 day. 

For fluid-fuel reactors the problems of burnup 
are solvable by straightforward engineering. 

For fixed-fuel reactors the rapid attrition of the 
fuel profoundly affects reactor design, both (1) be- 
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cause it is desirable to have operating time long 
compared to refueling time and (2) because of the 
changes in core characteristics due to having a time 
variation in fuel and poison concentrations. Conven- 
tional techniques [as used at the Materials Testing 
Reactor (MTR), the Engineering Test Reactor (ETR), 
and CP-5, for example| embody a refueling operation, 
performed after the reactor has been shut down and 
permitted to cool off, in which the fuel elements are 
replaced gradually and one at a time. The reloading 
shutdown takes 2 days to a week, depending on the 
cooling time adopted, the mechanisms available, and 
safety restrictions adopted. Thus cycle times of 20 
to 70 days are necessary to achieve 90% operating 
time. 

Even with a short time —20 days—and a 100-Mw 
power, 2 kg of fuel is burned. Because thermal re- 
actors are characterized by critical masses of 1 to 
5 kg, this means that something of the order of an 
extra critical mass must be provided initially. (Ac- 
tually, it is desirable economically to provide ini- 
tially two or three extra critical masses so that core 
fabrication can be amortized over a longer period.) 


Control System. The large overload of fuel re- 
quires that there be an elaborate control system. The 
fresh core must be held down, either by burnable or 
movable poison or by controlling leakage. Since leak- 
age neutrons are what the experimenter wants, the 
latter method should be preferred, especially since 
core poisons additionally restrict neutron availability 
per unit power. Unfortunately, reactor designers have 
so far chosen the “conventional” approach of con- 
trol rods rather than the “unconventional” one of 
reflector motion—to the experimenter’s disadvan- 
tage. 

To fulfill the requirements of core life, several 
expedients are, fortunately, available. In addition to 
the use of fluid fuel, which has several obvious ad- 
vantages for the long run, two points may be capi- 
talized: 

1. Intermediate- and, particularly, fast-reactor 
cores show relatively less reactivity swing over a 
large fractional core burnup. This not only makes 
control generally easier but also makes it more 
feasible to adopt leakage control. Moreover, such 
reactors do not show large xenon-buildup effects 
when shut down briefly. 

2. There is no physical reason for not having 
brief shutdowns for partial refueling, or, in the limit, 
designing for on-line refueling. The required engi- 
neering development to permit this, and to satisfy 
appropriate safety requirements, must ultimately be 
done. 


CRITICALITY AND POWER SHAPE 


Up to this point we have assumed that the high-flux 
reactor core can be made critical with an assumed 
volume, reflector material, configuration, and low 
parasitic core absorption. But this is not always the 
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case. Nearly all reactor types of potential utility can 
be made critical in volumes of 100 liters or more. 
However, attempts to capitalize on geometric factors 
by going to smaller cores require a progressive 
restriction of reactor type, and the smallest cores 
must be either proton-moderated or fast —a fact well 
illustrated in the choice of concepts for development 
as space nuclear power reactors. 

Small cores pose another problem. They tend to 
be relatively black to low-energy neutrons. There- 
fore, when immersed in a pure moderator, they see 
a large reentrant low-energy flux, which creates a 
power spike at the interface. Although this is desir- 
able experimentally —it produces neutrons close to 
the experiment —it is undesirable from the viewpoint 
of the core designer on three counts: (1) the power 
spike is a region of unusually high burnup and thus 
is not constant over life, (2) it limits total system 
power because the maximum power density is at the 
spike, and (3) the spike is likely to consist largely of 
epithermal fissions, the energy range least produc- 
tive of experimental neutrons. 


Consider New Concepts 


The most compact reactor systems that can be 
made critical are the proton-moderated thermal re- 
actors and fast reactors. In view of our familiarity 
through experience with them, water-moderated and 
-cooled reactors should be considered first. 


SMALL-CORE H,0 REACTORS 


There exists a body of information on the critical 
properties of these systems with H,O reflectors 
(p. 182 ff in Ref. 7). As practical solutions, with fis- 
sion products, *“*U, etc., water-reflected reactors 
can be made to go critical ina variety of sizes. Cores 
from 20 to 100 liters in volume all would require < 3 
kg of **U and would be quite thermal in their power 
spectra. A small core might have a 20-liter volume, 
ak, of 1.68, and a *U concentration of 100 g/liter; 
a middle-sized one would have V = 50 liters, k,, = 
1.40, and c = 36 g/liter; a large one, V = 100 liters, 
Rew = 1.25, and c = 25 g/liter. 

These cores would all be power-density limited. 
The “figure of merit,’’ power times availability fac- 
tor times geometry factor, could be approximated by 
[volume/(7 + 6)*] [(k,,-1)/k,,]. The results for the 
three sizes are almost identical. 


Beryllium Reflector. The substitution of beryllium 
for H,O as reflector would reduce critical-core radii 
by 6 to 7 cm. The result would be to permit a critical 
volume of ~20 liters at a k. of ~1.39 and a volume 
of ~7 liters at k., of ~1.7. 

The smaller beryllium-reflected system may be 
profitably examined. For one thing, high power densi- 
ties are more readily attained in small sizes than in 
large ones. Moreover, more concentrated systems 
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can tolerate structural neutron poisons more readily 
than can dilute ones. 

Although I have computed 7 liters as critical vol- 
ume, reduced water density at power makes 10 liters 
a safer volume for a reactor of k, = 1.7 and c = 100 
g/liter. At 10 Mw/liter and a flow area of 50 cm? in 
the outlet pipe, this would have, as a homogeneous 
reactor with allowed temperature rise of 200°C, a 
coolant flow velocity of ~25 m/sec in the piping. The 
realization of such a reactor would probably be as 
an aqueous slurry. 

Fluxes at beam locations would be several times 
10%. The major problems would be slurry circula- 
tion, radiolysis, and radiation damage to the con- 
tainer. It is probable that these could only be solved 
by using disposable (or, at least, dismountable) quick- 
connect primary components. 


HEAVY-WATER REACTORS 


A crucial difficulty with light-water reactors is 
that, unless they are quite small, the critical mass 
of fuel must compete with a significant quantity of 
water absorber for neutrons. Systems of more re- 
laxea size—100 liters or more —have a low Rk. and 
thus a small availability of excess neutrons per unit 
power. If one wants to keep the bulk of the fissions 
in the thermal range and yet have a large-volume 
core, the obvious step is to replace light water by 
very weakly absorbing heavy water. 

A practical heavy-water core, with some fission 
products and xenon, could go critical at a 100-liter 
volume with ak, of 1.88 and a ?*U concentration of 
about 6 g/liter. This is as estimated for a D,O re- 
flector, which is a convenient experimental medium. 
If, in some way, 1000 Mw could be extracted, the 
mean thermal flux would be greater than 4 x 10%, 
and the peak would approach 10”, 

The core would, as a sphere, have a radius of 
slightly under 30 cm. At an experimental point 40 cm 
from the core center, the current would be about 
1.8 x 10" neutrons/(cm?)(sec), and the small-beam 
flux would be close to 6 x 10%8, 

The great improvement over the light-water sys- 


tem comes from several factors: (1) the artificial 
one of this article —had I chosen 200 Mw at 10 Mw/ 
liter or 100 Mw at 5 Mw/liter, 20 liters would be the 
optimum-size core (although D,O systems can be 
made critical at this size, their concentration ap- 
proaches that of H,O systems), (2) the higher ,, of 
D,O cores makes them more productive of experi- 
mental neutrons for a given power, and (3) the higher 
flux/current ratio of D,O reflectors, which again be- 
comes apparent only in large sizes. 

In summary to this point: D,O systems with core 
sizes up to 100 liters have intrinsic neutronic advan- 
tages over H,O reactors. If 10 Mw/liter core power 
can be attained, D,O reactors are very attractive. 

The low critical concentration—cited as 6 g/ 
liter —makes it feasible to reconsider solution cores 
for this concept. 








158 POWER REACTOR TECHNOLOGY 


Central Zone. The very high core fluxes also make 
it attractive to consider locating experiments in an 
unfueled central zone. The peak core thermal flux 
for a constant-power reactor occurs with the core of 
minimum critical mass. According to my estimate 
this would occur with a core of ~120 liters—only 
slightly larger than the one just postulated —at which 
size the critical concentration would be 5 g/liter. A 
10" flux would be attainable in a central D,O tube in 
such a reactor. (Not yet a flux trap; in D,O it takes 
a very “gray” annular core and a large gap to develop 
maximum flux.) 

The flux in a homogeneous reactor can of course 
be improved further by on-site fuel processing to 
remove fission products, with a consequent decrease 
in critical mass; this can be done for this reactor. 
There is only one other engineering problem that 
was not mentioned in the discussion of high power 
density from homogeneous reactors—the problem 
that all structural parts, including core pot and en- 
trance pipes, must be very weak neutron absorbers. 
Otherwise, the virtues of the low absorption in D,O 
would be seriously compromised. Materials such as 
beryllium, BeO, carbon, and magnesium, with zir- 
conium as a possible cladding material, must be used. 

As a possible alleviation of the high pipe flow rates, 
an “octopus” or “dodecapus,” as sketched in Fig. 3, 








= 


cd 


¥) 


Fig. 3 Octopus-style homogeneous reactor would distribute 
out-of-core flow over four loops. 
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might be used. Such a reactor would have neutron 
production to some extent in the arms; and I can 
think of no easy way to figure out its real critical 
properties. 


SMALL FAST REACTORS 


Concentrated fast reactors are the smallest that 
can be made. For example, the critical mass ofa 
2351) sphere in a tub of H,O is only about 23 ke (p. 582 
of Ref. 7), corresponding to a core volume of 1.225 
liters and a core radius of 6.64 cm. In D,O the criti- 
cal radius is ~5.6 cm. Even a bare sphere (such as 
Godiva) has a radius of only 8.4 cm. Plutonium 
spheres are even smaller; a high-density plutonium 
sphere in a thick beryllium reflector has a critical 
radius of < 3.2 cm. 

These systems thus approach the “point source,” 
and so the possibility of capitalizing on geometry 
must again be examined. 

Selecting a point at radius 15 cm in a beryllium 
reflector as our experimental location, and assuming 
that a plutonium source has vy = 2.9 and k.=2 (to 
account for the rather large fraction of epithermal 
fission), we have a net neutron availability of 4.7 x 
10'* neutrons/Mw, a neutron current of 1.66 « 10'°/ 
Mw at the experiment, and a thermal-neutron flux of 
almost 3 x 10'4/Mw. 

A conservative example of the virtue of this sys- 
tem is illustrated below. It is conservative because 
performance could be markedly improved through 
use of plutonium and incorporation of very primitive 
internal cooling. 

Let us consider a concentrated uranium sphere 
immersed in a 5-cm-thick sphere of H,O. The criti- 
cal radius would be under 7 cm, and the total as- 
sembly radius’ would be 12 cm. The thermal con- 
ductivity of uranium is good enough that this assembly 
could operate at 8 kw for a uranium temperature rise 
of 200°C. At this power, given k, = 2.2, the beam 
current from the water surface would be 1.9 x 10'! 
neutrons/(cm2)(sec). In a graphite hohlraum this cur- 
rent may be somewhat intensified. This beam current 
is about what one would get from a 10” flux research 
reactor —admittedly, not a high-flux system, but ex- 
tremely good performance when compared with a 100- 
kw swimming pool! When such a reactor is opened 
for cooling, it can still be quite small. We con- 
sider now the next largest size, that of the Experi- 
mental Breeder Reactor I (EBR-I)® and Clementine.’ 

The EBR-I, as an example, has a core of a little 
over 6 liters, corresponding to a sphere of ~11-cm 
radius. With a most conservative cooling system, it 
produced 1 Mw of power. Its k,, was about 2.1. 


Sodium Coolant. We now know how to achieve power 
densities of 1 Mw/liter in large sodium-cooled fast 
reactors, and so it would be straightforward engi- 
neering to do the same for a 6-liter fast reactor. The 
use of beryllium instead of uranium for the blanket 
would markedly increase reactivity, making it pos- 
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sible to use smaller rods and more sodium in the 
core. So a power density higher than 1 Mw/liter can 
be predicted for this type of system. Sodium coolant 
flow of 15 m/sec seems possible, and, at 30% core 
porosity and 10-cm core radius, this would amount 
to a mass flow rate of 115 kg/sec. Because heat 
capacity is ~1.25 joules/(g)(°C), a temperature rise 
of 400°C would permit a power of 57.5 Mw, or almost 
10 Mw/liter. 


Short Core-Residence Time. The key to this high 
power density is the small core size, which permits 
the coolant to have a very short core-residence time. 
Of course, all the numbers cited are advances over 
current practice (EBR-II, for example, has a core 
2.5 times as long, sodium velocity % as great, and 
temperature rise A to he as much, in achieving some- 
what less than 1 Mw/liter). Nevertheless, the sodium 
route is one of the more straightforward ways of 
achieving high power density. 

At 60 Mw in a 10-liter *U-fueled core witha 
beryllium reflector, we might achieve k.. = 2.1. With 
an experimental position at 24 cm, the experimental 
current would be 3.4 x 10", and the small-beam flux, 
around 7 x 10%, 


Hydrogen Coolant. A somewhat larger core could 
undoubtedly be built of refractory materials. We have 
investigated tungsten—UO, cermet fast-reactor sys- 
tems for nuclear-rocket application. At about 50 
liters and up, it becomes possible to extract 1000 
Mw from them when the coolant is hydrogen, which 
can achieve high flow rates in small passages with- 
out unusual pressure drop and has a very high heat 
capacity per unit weight. The refractory core per- 
mits the heat source to sustain a high driving tem- 
perature difference for heat transfer to the bulk 
coolant. 

Such a system — 1000 Mw, 20-cm radius, k.. ~ 1.8 — 
would yield the impressive flux of ~ 1.6 x 10". Un- 
fortunately, it is probably too costly to supply many 
cores of this type, and the cost of using hydrogen as 
a coolant is also extreme —as is likewise the hazard. 
Nevertheless, it indicates that technology in itself 
does not limit power density. 


CAVITY AND THIN-SHELL REACTORS 


It has been often observed, in this article and else- 
where, that a major limit to power density in hetero- 
geneous reactors is a relatively long coolant-flow 
path through the core. Pressure drop increases with 
length, and temperature rise is constrained to some 
maximum. It is for these reasons that smaller cores 
may have higher power densities. 

It is clear that shorter coolant path lengths can be 
achieved with radial coolant flow. For a more-or- 
less compact core, only a factor of 2 reduction in 
flow length is possible, and this is hardly worth the 
design difficulty. However, there are core types for 
which the fuel configuration can be a thin annular 


section. For such cores, quite short flow paths can 
be realized by radial flow. 

The most intriguing of these systems are the D,O 
flux-trap cores and their close relatives, D,O cavity 
reactors. Several reactor systems of this general 
sort were investigated some years ago.'° The basic 
observation is that, in a large body of moderator 
(D,O being the best example), a spherical shell of 
fuel just thick enough to be black to thermal neutrons 
will be critical.‘ If the shell is made larger, its 
blackness can be reduced to some limiting value, 
which is that required for criticality in a thin slab. 
Moreover, if the interior of the shell is evacuated, 
only a slight increase in blackness keeps the system 
critical; and, if the contents of the shell are dis- 
persed uniformly within the sphere (the cavity reac- 
tor), again only a little more fuel is required. 


Cavity Reactor. The cavity reactor is often con- 
ceived in terms of a fissioning gas with temperatures 
approaching those of plasmas. It is a possible core 
for plasma-reactor schemes, such as those of Refs. 
11 to 13. The cavity reactor may also be fueled with 
a suspension of fissioning “dust”; as such, it is ef- 
fectively a gas-phase homogeneous reactor of poten- 
tial quite similar to that of the large homogeneous 
heavy-water solution. Problems of corrosion and 
plateout are replaced by problems of erosion and 
adhesion. 


Thin-Shell Reactor. The existence of thin-shell 
critical assemblies presents a more novel reactor 
concept. Figure 4 illustrates a possible configura- 
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Fig. 4 Colander type D20 thin-shell reactor. 
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tion: The core is a thin shell of a strong fueled ma- 
terial—such as a stainless-steel—urania cermet. 
This shell is perforated with a very large number of 
small holes, possibly drilled or punched at the time 
of fabrication, and probably clad by a chemical or 
electrochemical technique. The shell has at least 
one reentrant thimble for high-flux irradiations and a 
stem for coolant entry. Coolant is pumped in, forced 
through the holes of the colander, and dumped into 
the main moderator vessel . 

Although the inlet pipe must contain the high flows 
characteristic of homogeneous reactors and the shell 
must support significant pressure, the high-tempera- 
ture outlets of the homogeneous reactor have disap- 
peared, and pressure drops do not approach those of 
parallel-flow heterogeneous systems. The D,O may 
be permitted to evaporate and superheat in passing — 
it will recondense in the outer tank. (Not shown is 
the main cooling system, which withdraws warm 
moderator from the tank and returns part of the cold 
moderator by downward circulation around the out- 
side of the shell.) 

For those who are concerned with the strength of 
such a colander, an alternate scheme has been pro- 
posed.‘ The reactor shell is now solid (Fig. 5). The 
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Fig. 5 Sprinkler type D20 thin-shell reactor. 
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coolant is dispersed by flow through a nozzle head 
and impinges on the inner surface of the shell. The 
velocity of the coolant jets is high enough so that they 
can cut through the steam which comes off that sur- 
face, impinge on it, and be evaporated. The outer 
surface of the shell would be vapor bound in stable 
film boiling but not at a high enough temperature to 
burn out. 

Finally, the scheme may be reduced to more con- 
ventional terms by reverting to parallel flow in small 
tubes of some appreciable length. Such a reactor 
looks like the colander of Fig. 4, with the thin shell 
replaced by a tubular fueled structure of 2 to 10 cm 
thickness (and flow length). 

The thin shell may also be used in reactors of 
more normal flow design. D,O flux traps may be 
simulated by heterogeneous designs consisting of 
individual fuel elements spaced around a cylindrical 
annulus or a pseudospherical “basket.” The large 
radius of optimal D,O flux traps permits, in fact, 
enough space for a pressure-tube reactor design, 
which might make on-line reloading more feasible. 

Another version of a thin-shell reactor is one com- 
puted by Ergen.*»® It seems possible to build a beryl- 
lium-moderated version in which the shell contains 
circulating molten salt. The excellent characteristics 
of molten salts as possible homogeneous reactor 
fuels (high permissible temperature rise and high 
gravimetric heat capacity) make a 100-Mw power for 
such a system credible. Ergen’s calculation predicts 
a7 x 10% central-island flux for this case. 


This article is based ona paper presented at the 
Symposium on Intense Neutron Sources,!® Santa Fe, 
N. Mex., Sept. 19-23, 1966. 
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PHYSICS 


Zone-Loaded Critical Experiments 
for Fast Breeders 


By George J. Fischer, Armando Travelli, and Daniel A. Meneley 


Physics uncertainties—-among others—limit the 
technology of the large fast power breeder, a reactor 
that might be fueled with 2000 kg of plutonium ina 
6000-liter core and be capable of 1000 Mw(e). Some 
of these uncertainties influence safety analyses. 
Others relate to the behavior of the nuclear fuel 
over the fuel cycle and influence calculations of 
such parameters as doubling time and reactivity 
Swing over core lifetime. These uncertainties force 
appropriate conservative design decisions that tend 
to lower the economic potential of the fast breeder. 

Experiments on zero-power critical assemblies 
can effectively test the degree of uncertainty of 
many important calculations. They can also provide 
useful basic data to help reduce the uncertainties. 
There are some difficulties involved in straight- 
forward mockup experiments, however. These dif- 
ficulties will be discussed here; it is largely because 
of them that zone-loaded experiments are important. 

A cross section of the material arrangement of a 
typical zone-loaded critical experiment is shown in 
Fig. 1. It is a zero-power critical assembly! and 
contains only a small fraction of the full core of 
interest. This fraction of the full 7**Pu—***u fueled 
large breeder core is mounted as a center zone 
inside a “driver” region that supplies the additional 
reactivity needed to achieve criticality. A “buffer” 
or “filter” zone may be used between these two 
regions to improve such characteristics as the 
neutron-energy spectrum near the edge of the center 
zone. The driver is fueled with 7°U, which is more 
readily available than plutonium. 

This article presents the reasons why zone-loaded 
experiments are important and describes a number 
of physics questions that can be answered through 
these experiments. It also discusses the limitations 
of the zone-loaded method, the safety aspects of this 
type of critical loading, and the major needs for 
further development of the concept. An example is 
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Fig. 1 Cross-section view of zone-loading assembly. Re- 
actor materials are often loaded into matrix tubes of ~2- 
by 2-in. cross section, which gives the rectangular pro- 
file.’ 


given of a critical test of the ability of the analyst 
to predict full-core reactor properties on the basis 
of analysis of a zone-loaded critical experiment. 


Why Zone-Loaded Experiments? 


Among the important advantages of zone-loaded 
experiments are: 

e Smaller and presently available reactors such 
as ZPR-VI (Fig. 2) can be used effectively. 

e Much smaller inventories of materials such as 
39pu or 73y fuels are needed in many cases. 
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e Cost of items directly related to reactor size is 
reduced. 

e A program of survey type experiments on zone- 
loaded reactors could provide important basic data 
for a range of breeder-reactor compositions, possibly 
years sooner than these data might otherwise be 
obtained. 

e There may be Safety and health advantages to 
the use of zone-loaded mockups. 


Difficulties involved in straightforward full-core 
mockups are numerous. The 6000-liter core, for 
example, would require a very large investment in 
fuel plates and other mockup materials —in amounts 
that usually are not obtainable at present. (Many 
nations will be forced into doing zone-loaded experi- 
ments because of limited availability of these ma- 
terials.) Even when the materials are on hand, from 
4 to 8 weeks would be needed just to assemble one 
full-core critical experiment. The radioactivity from 
*40Du, fission products, and possibly ***U and its 
associated isotopes, all or some of which might need 
to be included, would be a hazard to operating person- 
nel. Decay heat from these materials would induce 


reactor drift, which could limit the precision of 
reactivity measurements. Because of the high operat- 
ing costs for facilities in which zero-power critical 
assemblies are built and because of the pressing 
need for experimental data, especially sodium-void 
and Doppler-effect measurements, it is essential 
that these facilities be employed with maximum 
effectiveness. 


Physics Needs of the Fast Power-Breeder 
Reactor Program 


Phenomena of particular concern in the difficult 
area Of safety of the sodium-cooled large fast 
breeder—and comparatively difficult to calculate 
accurately —are the sodium-void effect and the Dop- 
pler effect. For some relatively simple and typical 
accident calculations that include both effects, sodium 
expulsion has been calculated to add reactivity at 
rates of 2 to 100 dollars per second, Addition of 
1 dollar of reactivity to a reactor operating ata 
steady state (delayed critical) makes the reactor 
prompt critical—a dangerous condition in which the 
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power doubles every few thousandths of a second. It 
is very important that careful experiments be per- 
formed to check the ability of the theoretician to 
calculate these two phenomena. In addition, since 
these effects have been calculated to be highly 
sensitive to the presence of higher plutonium iso- 
topes, fission products, and the nature of the clad 
materials, the experiments must include these ma- 
terials. 


THE DOPPLER EFFECT 


The Doppler effect is aprompt shutdown mechanism 
that operates primarily by absorbing neutrons in the 
energy range from about 1 to 15 kev, the low-energy 
tail of the neutron-energy spectrum of a fast reac- 
tor. This absorption takes place primarily in the 
fertile atoms, which exceed the fissile atoms by 
ratios of about 7:1 to 9:1, in the fissile isotopes, 
and, to a smaller extent, in structural materials. 
The effect occurs as a response to temperature 
changes, as in the heating of the fuel during a power 
surge. An objective of fast reactor designs is to 
achieve a Doppler-effect shutdown mechanism suf- 
ficiently strong that possible gains of positive re- 
activity upon voiding of sodium would be safely 
controlled. 


THE SODIUM-VOID PROBLEM 


Sodium occupies 50 to 70% of the core volume in 
the more common fast power-reactor designs. A 
relatively effective moderator because of its modest 
atomic weight, sodium shifts the neutron spectrum 
toward lower energies. In addition, sodium also has 
a large and broad resonance that peaks at 2.85 kev 
and strongly affects the neutron population in the 
1- to 4-kev energy range, which contributes to the 
Doppler effect. Full or partial loss of sodium cool- 
ant, such as might follow pump failure or a transient, 
reverses these spectral effects. There is thus a 
spectral component of the reactivity change associ- 
ated with the loss of sodium. 

This spectral change which accompanies sodium 
voiding usually adds reactivity. Loss of sodium ina 
region about the reactor core center of some large 
fast power-breeder designs can increase reactivity 
by several dollars. On the other hand, loss of sodium 
from any reactor region in which there is a spatial 
flux gradient (as near the core boundary) increases 
the leakage of neutrons from the core and thereby 
tends to reduce the reactivity in competition with the 
spectral component. This reactivity change is the 
other component of the sodium-void effect and is 
called the leakage component. 


OTHER NEEDS 


In addition to the sodium-void and Doppler effects, 
there are many other uncertainties to be resolved. 
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The theoretical physicist would like to check his 
cross-section choices for the materials in the reac- 
tor by means of material-worth measurements and 
of neutron-energy-spectrum measurements. He would 
like to check the reliability of some ofhis simplifying 
assumptions and procedures, such as the spatial 
dependence of his cross sections, heterogeneity cor- 
rections, choice of transport approximations, coupling 
between modules for a modular concept, and the 
change of this coupling with reactor operation. There 
are engineering problems that are sometimes dif- 
ficult to calculate or important to verify. These in- 
clude worths of partially inserted rods and effective- 
ness of shielding, for which additional approximations 
must be introduced and checked. 

Plutonium became available for zero-power fast 
critical experimentation only a few years ago. Be- 
cause physicists associated with the critical as- 
semblies employing plutonium were required to de- 
vote a large part of the assembly availability time 
to mockup experiments, very little good basic ex- 
perimental work on the reactor physics of plutonium — 
uranium fueled fast reactors exists. 


Usefulness of Zone-Loaded Experiments 


The central zone of the assembly! shown in Fig. 1, 
Assembly 45A of ZPR-III, mocked up a large 7°*Pu— 
387 monocarbide fast breeder reactor very satis- 
factorily with the use of only 51 kg of 7°Pu. The 
21-cm radius of this central zone may be close to 
the lower limit of what is desirable on a zone-loaded 
experiment. A cylindrical geometry was used, with 
an axial core length of 91 cm chosenas a representa- 
tive fast breeder core height. The zone-loaded as- 
semblies discussed here will be assumed to have the 
arrangement of this assembly. 


MEASURING NEUTRON-ENERGY SPECTRA 


Many physics measurements related to the neutron- 
energy spectrum of large fast breeder reactors can 
be performed satisfactorily in a zone-loaded reactor. 
The most direct measurement is that of central neu- 
tron spectra. The experiments performed by Ben- 
nett, Gold, and Huber’ provide an excellent example 
of the accuracy and resolution of the spectral data 
that can be obtained from a zone loading. Valuable 
energy-spectrum measurements in the direction of 
the axial blanket are equally possible in zone-loaded 
cores, 


MEASURING SODIUM-VOID EFFECTS 


Sodium-void measurements of several types can 
be performed satisfactorily on zone-loaded reac- 
tors,'°*4 The spectral component of the sodium -void 
effect is easily measured at the center of the zone 
loading, as are variations in the spectral effect due 


Fall 1966 PHYSICS 165 


to composition changes. Also straightforward are 
measurements of the important effects of channeled 
loss of sodium (related to the accidental voiding of a 
fuel subassembly). The leakage component of the 
sodium-void effect in the axial direction may also 
be measured. These experiments should be limited 
to one or more channels around the axis of the re- 
actor. They frequently are sufficient to verify theo- 
retical methods. 

Care must be taken, of course, that the flux 
perturbations due to any measurement, especially if 
made at some radial distance away from the axis of 
the central zone, do not extend to a region where the 
neutron spectrum is affected by the buffer and driver 
zones. Although these effects might be compensated 
for to a considerable degree, it is always desirable 
to keep the measurement as free from uncertainties 
as possible. 

As another precaution, the analyst must remain 
aware of the possibility that, in the approach of 
first-order perturbation theory, the denominator term 


f 
f (= ohn) x (= 02,9) dVol 
j ; 
may change upon removal of sodium over a significant 
region of the central zone. 


MEASURING THE DOPPLER EFFECT 


Doppler-effect measurements have been made at 
or near the center of the zone loadings of the ZPR- 
III reactor on Assemblies 43, 43A (Ref. 5), 45, and 
45A (Ref. 1) and in zone-loaded Assembly 4Z of 
ZPR-VI at Argonne.’ Atomics International has per- 
formed Doppler-effect measurements in zone-loaded 
assemblies that employed thermal drivers.’ All these 
measurements used small samples, which allowed 
much freedom in such matters as the nature of ma- 
terials that surrounded the samples to be heated. 

Experimental difficulties in Doppler-effect mea- 
surements made thus far in zone-loaded assemblies 
are only in small part different from the difficulties 
that would be found in equivalent measurements on 
nonzoned assemblies. Particular care must be taken, 
however, in attempting to measure the change in 
Doppler reactivity effects when sodium is removed 
from the environment surrounding the heated sample. 
It is very easy for the domain in which the sodium 
is removed to become sufficiently large that the 
Doppler element “sees” the presence of the buffer 
zone. 

With moderate care, multiple-sample and “pot” 
type Doppler-effect measurements should be suc- 
cessful in zone-loaded assemblies. 


KINETICS STUDIES 


Rossi-alpha type kinetics measurements can be 
performed in zone loadings and apparently can be 


interpreted to predict the full-core kinetic character- 
istics with considerable confidence.® This property 
was demonstrated first in ZPR-III experiments on a 
coupled fast thermal-reactor concept,’ a form of 
zone loading. 


EQUIVALENCE TO MODULAR FAST 
BREEDER SYSTEMS 


The physics of the modular concepts can be related 
clearly to the physics of the zone-loaded concepts. 
The radius of the central zone can be increased so 
that its reactivity grows, and the nature of the buffer 
zone can be changed to approximate the blanket 
region between modules. By letting the driver region, 
now with perhaps the same fuel as the central zone, 
represent a simplification of the modules external 
to the central module, we achieve a “zone-loading” 
representation of the modular concept which is 
particularly amenable to analysis. Much ofthe statics 
and dynamics of the physics of interaction of the 
modules could be checked by this type of experiment. 
For instance, the change of the coupling between the 
modules with coolant loss or fissile isotope buildup 
in the blanket could be measured readily. 


EQUIVALENCE TO FUEL-IRRADIATION FACILITIES 


Some proposed fast reactor fuel-irradiation facili- 
ties are, in a sense, zone-loaded reactors. They 
frequently employ a central region containing fuel 
elements whose irradiation characteristics are to be 
learned, a buffer region to shape power distribution 
or neutron-energy spectrum, and a driver fuel chosen 
largely for its power density, burnup, and fabrication 
characteristics. Thus techniques learned during zone- 
loaded critical-assembly design can be useful to 
designers of fuel-irradiation facilities. 


Limitations of Zone-Loading Procedures 


A not unimportant difficulty of the zone-loading 
procedure as compared to the single-zone procedure 
is that the physicist must be more careful about the 
design and interpretation of some experiments. For 
example, the removal of sodium from a large portion 
of the central zone is an experimental variation that 
must be attempted only with caution. For such an 
experiment the best choice of buffer-zone composi- 
tion becomes a matter of concern. For certain 
regions, transport-theory calculations may predict 
spectrum distributions significantly different from 
those predicted by diffusion-theory calculations. 
Greater care must be taken in the procedure for 
collapsing multigroup cross sections from ultrafine 
data because the spectrum over some of the groups 
may change considerably within short distances. 
However, awareness of these difficulties should as- 
sure satisfactory interpretation of zone-loaded ex- 
periments. 
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ESTIMATING REACTIVITY CHANGES 

Another difficulty arises in trying to predict from 
the result of a reactivity measurement in a zone- 
loaded reactor the result that the same measurement 
would have yielded if it had been made in the full 
reactor that the center zone of the zone-loaded re- 
actor represents. A simple example of this difficulty 
is found when a reactivity measurement is to be 
made at the center of the zone loading and when a 
perfect match of real and adjoint neutron-energy 
spectra with the corresponding spectra in the single- 
zone reactor has been achieved over the region of 
the reactor affected by the perturbation. A satis- 
factory approximation to this condition is achieved 
relatively easily with careful design of the zone 
loading. 

If this measurement were to be interpreted by 
means of first-order perturbation theory, the ex- 
pression for the perturbation of k would be 
Ak = el o* Alx,v, Zi, fe Togs a One | p, dVol 
z Roasts gt [x,y Dy; /k| >, dVol 
where all terms have their usual definitions. Both 
integrals extend over the entire reactor, but, of 
course, the terms in the numerator are nonzero only 
over the regions in which the perturbation is made. 
With the assumption that the real and adjoint fluxes 
in the zone loading match those of the full core over 
the region of the perturbation, the only difficulty in 
relating worth measurements made at the center of 
the two cores lies in the calculation of the denomi- 
nator terms for the two reactors. For the full core 
this calculation is relatively straightforward. For 
the zone loading, each of the zones adds several 
sources of uncertainty. The results of an investiga- 
tion of some of the more difficult calculational 
aspects of this problem are given in the subsection 
entitled “A Test of the Zone-Loading Method.” 


ESTIMATING CRITICAL MASS 


The critical mass of a large fast power breeder 
cannot be directly determined by a zone-loading ex- 
periment. However, this is not necessarily a severe 
limitation. Many large power breeders will employ 
several zones of different enrichments to provide 
power flattening. Small changes in the location of 
the interfaces between the zones will result in 
criticality of the actual reactor. Many important 
indirect questions (e.g., material worth and detailed 
neutron-energy spectra at the core center and near 
the axial blanket) used to check the ability to cal- 
culate critical mass can be answered very well by a 
zone loading. Again, recent cross-section measure- 
ments and refinements in multigroup cross-section 
preparation procedures have improved the ability to 
estimate critical mass. Whereas errors inpredicting 
critical mass were 30 to 40% for difficult cases a 
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year or two ago, errors larger than 10% may be 
expected to be relatively unusual in the near future. 
The problem of predicting the critical mass of a 
large reactor by zone-loading procedures has been 
treated with some success recently in France." 


Safety Aspects of Zone-Loaded 
Critical Experiments 


It would be thoroughly misleading to pretend that 
the issue of safety of zone-loaded critical experi- 
ments, even relative to that of full-core loaded re- 
actors, can be adequately treated in a few sentences. 
The statements made here are intended only to point 
out important aspects of a complex problem. 

A full-core (nonzone loaded) loading on a zero- 
power critical reactor can require from several 
hundred to several thousand kilograms of highly 
toxic plutonium fuel. If, in spite of the very con- 
siderable precautions taken in the design and opera- 
tion of these reactors, a severe accident were to 
occur (definitely a highly improbable event), the 
containment of the reactor might be ruptured and 
very dangerous materials dispersed for several 
miles downwind. For this reason, operators of fast 
reactor critical experiments give much attention to 
safety. 

An obvious aspect of the relative safety of the two 
types of experiments is that the zone-loaded reactor 
would contain 50 to 200 kg of ?**Pu in the central 
zone, approximately one-tenth as much as the full- 
core reactor mockup. The driver region of the zone- 
loaded core usually would be fueled with 200 to 
300 kg of the far less toxic ?®U, 


SHUTDOWN MECHANISMS 


An important prompt shutdown mechanism of the 
critical assembly comes from the Doppler effect. 
This is achieved if the fertile °°U is mixed with the 
*39Pu in a ratio of about 3:1 when the fuel plates or 
rods are fabricated and if the assembly is composed 
of such other materials that the real and adjoint 
neutron-energy spectrum in the kiloelectron-volt 
region is sufficiently large. This second requirement 
may not be satisfied. It is especially difficult to 
provide an adequately soft spectrum in the driver 
region of a zone-loaded assembly. Furthermore, the 
driver fuel of older reactors may consist of iy 
separated from the 7**U, Under these conditions the 
23517 will undergo many more fissions than the se | 
and therefore will heat more rapidly during an inci- 
dent. Any Doppler effect in a driver of this type 
might add to the reactivity rather than subtract from 
it because of the generally positive nature of the omey 
Doppler effect. 

Separated driver fuel, however, has better re- 
activity-loss characteristics through expansion on 
heating than does mixed driver fuel. 
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How well can reactor analysts use data from a 
zone-loaded reactor experiment to predict the re- 
sults for the same measurements in the full non- 
zoned core? We were recently given the opportunity, 
and challenge, to make such an analytical study and 
prediction. Although the success of the analytical 
prediction made for our test experiment does not 
assure the credibility of all analyses of zone-loading 
experiments, the study demonstrates that successful 
predictions of the properties of large fast breeder 
reactors can be based on zone-loaded experiments, 

We were asked to analyze zone-loaded Assembly 
4Z (in ZPR-VI) and to predict the results of the more 
crucial subsequent experiments on Assembly 5, a 
single-zone 2600-liter core of the same composition 
as the zone-loaded core. The assignment was re- 
solved into several parts. One study was that of the 
deviation from the ‘‘correct’’ answer which results 
from the simplifying assumptions an analyst fre- 
quently makes in the course of his work. Examples 
of these approximations are the use of diffusion 
theory vs. transport theory and the use of one- 
dimensional calculations vs. the much more costly 
multidimensional calculations. A second study con- 
centrated on the difficulties of evaluating the de- 
nominator terms of Eq. 1. Because of the nature of 
the construction of Assembly 4Z, it was evident that 
the accuracy of any prediction of Doppler, sodium- 
void, or other worth measurement would be sensitive 
to how well this term was calculated. The third ob- 
jective was to predict the critical mass of Assembly 
5. Much of the detail and physics of these studies 
has been published.!? Only a summary of the results 
is given here. 

A number of interesting results were obtained 
from the first study. It was shown that transport 
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theory predicted very different behavior of the neu- 
tron flux, in certain energy intervals, near the edge 
of the central zone and through the filter and driver 
regions. The 2%8U fission activity was significantly 
different, so that, if Assembly 4Z had been a fuel- 
irradiation facility, diffusion theory would have rather 
poorly predicted the power density due to 2°8U fis- 
sion. Other results showed cancellations of detail 
effects in Assembly 4Z, in some cases, which would 
help generate false confidence in the predictions of 
simpler methods. 


In the important area of predicting analytically the 
worths of measurements of reactivity made at the 
center of the full core from those at the center of the 
zoned Assembly 4Z, it was predicted that the ratio 
of these worth measurements would be a constant 
and that the value of this constant would be 1.003. 
(Heterogeneity effects in the driver were not ac- 
counted for in the prediction because of the com- 
plexity of the geometry used.) 


It was also predicted that the critical mass of the 
core of Assembly 5 would be 1664 kg of 2U, with an 
uncertainty of +150 kg mainly due to uncertainties 
in the cross sections used in the calculations. 


These predictions were tested when Assembly 5 
was built. It was found that the weighted average of 
the ratios of the reactivity worths of a large number 
of materials measured in the two assemblies was 
1.003 + 0.002, against the predicted value of 1.003. 
This agreement must be considered partly fortuitous. 
The experimental critical mass, corrected for edge 
effects, was close to 1580 kg of 2U, 84 kg less than 
predicted but well within the estimated uncertainty 
of the calculation. 








VIOLENT DISASSEMBLY 


The most improbable but also the most severe 
accidents are those for which all other accident- 
termination mechanisms have failed and the accident 
is shut down by violent disassembly. The accelera- 
tion of material is caused by gradients in the 
hydraulic pressure generated by heating of reactor 
material, especially the fuel, which may melt or 
vaporize. The local pressure in a region is usually 
proportional to the power density there. 


If a nonzoned large reactor is mocked up on the 
facility at the time of the accident, it may be the type 
that employs several regional fuel enrichments to 
obtain a flattened (economically desirable) power 
distribution. The flattened power distribution, how- 
ever, leads to very small pressure gradients in a 
severe accident. Thus disassembly and shutdown of 
nuclear multiplication processes would be achieved 
only after large amounts of nuclear energy had been 
released. 

A safety aspect of common zone-loading designs of 
the form represented by Fig. 1 is that they frequently 
are constructed, for practical reasons, with power 
densities in the driver zones which are from 1.5 to 
3 times larger than those of the central zone. Thus, 


during a severe nuclear incident, very large pres- 
sures would occur, most probably, in the driver 
zone well before they would occur at other locations. 
Also, there usually is a large gradient in power 
density across the driver in the direction away from 
the reactor center, a very desirable safety charac- 
teristic. We have studied the neutronic and hydro- 
dynamic behavior of several zone-loaded fast reactor 
experiments for severe accident conditions. The 
accident code Ax-1 (Ref. 11) and a more recent ANL 
code, PREAX,! able to handle meltdown histories as 
well as severe accidents were used. The calcula- 
tions show that initial high pressures in the driver 
region begin to drive the inner boundary of the driver 
toward the center of the core (a reactivity-adding 
mechanism) and at the same time to drive most of 
the remainder of the driver outward (a reactivity- 
losing mechanism). A very short interval of time 
over which the reactor gains a very small amount of 
reactivity was seen, followed by a relatively rapid 
shutdown. 

Because of the high power density of the driver, 
disassembly pressures in the driver region began 
early in the history of severe accidents calculated 
for zone-loaded reactors. Since small displacements 
of the driver produced rapid reactivity loss, nuclear- 
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energy release for accidents in zone-loaded reactors 
tended to be smaller than for comparable accidents 
in critical assemblies mounted with the equivalent 
full-core loadings. 

Broad generalizations of this type, of course, must 
not be used to conclude that all zone-loaded experi- 
ments are inherently safe. This problem deserves 
much more careful attention than it has received. 


Improvements Needed 


Limited inventories of nuclear materials, espe- 
cially isotopes of plutonium, make zone-loaded ex- 
periments the only available experimental means for 
many laboratories to study large fast power-breeder 
reactors. For others, the technique is an exception- 
ally useful research procedure. For these reasons, 
it is important to continue to extend the understanding 
and capabilities of the procedure. 

Some rather basic practical questions have not yet 
been adequately resolved. An example is the question 
of what minimal amount of plutonium fuel, containing 
a varied enrichment of higher isotopes, should be 
purchased to study satisfactorily the physics of large 
power breeders having this fuel composition. 


BETTER CROSS SECTIONS 


Particularly needed for successful analysis of 
zone-loaded assemblies is an analytical tool which 
would produce spatially dependent multigroup cross 
sections from ultrafine cross sections and which 
would take into account the flux-perturbing presence 
of materials in adjacent regions. One such procedure 
is being developed by D. Meneghetti, A. Rago, and 
K. Phillips at Argonne. 


HETEROGENEITY CORRECTIONS 

The coarse plate structures used as mockup ma- 
terials on zero-power critical experiments approxi- 
mate only roughly the structures of power reactors, 
and they create a problem that has been recognized 
for some years. A satisfactory correction is espe- 
cially needed for zone-loaded cores because driver 
regions are relatively thin and because they often 
include materials highly transparent to neutrons. 


EXTENSIONS AND TESTS 

There is a serious need to examine ways in which 
the zone-loading technique might be extended for 
even greater usefulness. There is also a need to 
identify clearly those cases in which zone-loading 
methods could lead to misinterpreted results. Finally, 
there is a need to specify a set of experiments that 
would most clearly test analytical procedures. 
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Further development of the zone-loaded reactor 
concept as an experimental research tool is im- 
portant to fast breeder-reactor development and is 
also an intriguingly complex problem for the analyst. 
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CONTROL AND DYNAMICS 


Let’s Apply Optimal Control Theory and 
Nonlinear Analysis 


By Ronald R. Mohler, University of New Mexico 


Reactor-control problems traditionally have been 
solved —or at least satisfied—by trial-and-error 
techniques with digital- and analog-computer simula- 
tions based on classical linear-systems theory. With 
the aid of recently developed | “modern”| systems the - 
ory, i.e., optimal control and nonlinear stability anal- 
ysis, we may achieve more precise solutions and gain 
the dual benefits of increased safety and economy. 

Rather than the theoretical details of optimal con- 
trol and nonlinear analysis, this article deals with re- 
sults and their implications, constructively criticizes 
certain areas, and indicates some of the pitfalls into 
which researchers have stumbled. 

Some of the easier reactor-control problems are 
discussed. You are most fortunate, though, if your re- 
actor fits such cases. Tough practical problems un- 
fortunately are the ones confronting reactor engineers. 
With the computations as involved as they are, it is 
no wonder the simplest possible problems are solved 
first. If the work is tobe appliedin everyday practice, 
however, systems people must dirty their sometimes 
sterile computations with practicality. Also, more 
cooperation is needed between the practical people 
and the theoretical people. Systems people must learn 
some practicality and reactor engineers must become 
knowledgeable in systems theory in order to break the 
communications barrier that often exists. 


Optimal Reactor Control 


The term optimum is by no means new to the prac - 
ticing reactor engineer. The designer is continuously 
trying to optimize such things as flux distribution and 
fuel loading. Indeed the practicing reactor-controls 
engineer often arrives at an intuitive “optimum” de- 
Sign by eyeballing a response while tweaking a pot 
setting or by examining various root loci. In the se- 
lection of such an optimum controller, various in- 
dexes of performance, such as response time, over- 
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shoot, damping, bandwidth and peaking, might be 
considered. 

The optimal control considered in this article, how- 
ever, is one that results from a more precise (al- 
though not necessarily more meaningful) mathematical 
problem. For example, the problem might be to vary 
reactivity in such a manner that reactor power is in- 


creased while some performance index, 
“ty 
J = L at 
24 


is held to a minimum. Here /, is an initial time, /; is 
a terminal time, and Lis the cost, which might be a 
function of power, temperature, reactivity, and other 
state or control variables. (State variables and state- 
space analysis are treated for the uninitiated reader 
in Refs. 1 and 2.) For certaincommercial reactors, L 
may be an economic cost in dollars per hour of oper - 
ation. In general, such problems have not been studied. 
The physical system further implies certain con- 
straints on the allowable optimal control. 

Optimal control problems have been solved by the 
maximum principle,‘ dynamic programming,°* ' and 
extensions of the classical calculus of variations.°”' 
The reader is referred to the indicated textbooks for 
excellent treatments of the theory along with ex- 
amples. 

Minimal-time control of nuclear reactors probably 
has been analyzed in more detail than any other opti- 
mal reactor -control problem? Time -optimal con- 
trol is a maximum allowable effort type of process (in 
the simplest form, a bang-bang control). Roughly this 
means that, to increase power level inminimum time, 
reactivity is added as rapidly as possible until some 
maximum allowable limit is reached—the whole 
process within values that guarantee safe operation. 
Since the control is not inertialess, reactivity must be 
decreased at some predetermined time so thatthe re- 
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quired reactivity is attained at the desired terminal 
power level. Further, precursor levels have not yet 
reached equilibrium, and the required terminal reac - 
tivity is time variant so as to maintain a constant 
power level. 

All this looks good on paper, but itdoes not work in 
practice. First, the model used to describe the sys- 
tem and used to arrive at the optimal process is al- 
ways erroneous. Second, the required terminal con- 
trol requires a precise measurement —estimation of 
precursor levels. Fortunately there are various ap- 
proximations to the time-optimal process which may 
be synthesized in practical ways. For example, when 
the appropriate time (or reactor power) is reached, a 
fast-correcting closed-loop process can be utilized to 
control the desired terminal value. Sucha closed-loop 
process can itself be designed optimally with a cost 
that is a quadratic function of power error and con- 
trol. Another effective technique utilizes a dither 
(multiple bang-bang) process with certain control rods 
jittering about their theoretically required positions. 
The switching is determined by a small allowable 
dead zone in reactor power .® 


STATE-SPACE CONSTRAINTS AND PRACTICAL 
SUBOPTIMAL CONTROL 

Complications in the optimal process also result 
from certain constraints in such state variables as 
power level (possibly reactor period) and tempera- 
ture. Although the computation of optimal processes 
with state-variable constraints is complicated by the 
presence of more Lagrange multipliers, it can be 
carried through with computers (as discussed in 
Refs. 8 and 12). Again the exact process is not too 
practical because of a lack of knowledge of the pre- 
cise system dynamics and state. As would be ex- 
pected, approximations can be used to synthesize a 
near optimal or suboptimal control system for the en- 
tire portion of the trajectory. A fast-acting closed- 
loop control can be utilized. Also, state-variable con- 
straints and control constraints can be handled by the 
proper power-demand limits (and any other control- 
loop inputs) as well as by the physical limits im- 
posed on control variables (e.g., rod-position demand 
limited). Another way to manage these constraints is 
to optimize the loop design with a performance index 
that implicitly includes constraints (e.g., as a certain 
quadratic performance index). For the nuclear rocket 
engine, a simple suboptimal control synthesized by 
cut-and-try methods with the aid of a computer simu - 
lation gives nearly optimal performance.*” A more 
exact optimal design is obtained in Ref. 13 by optimiz- 
ing the closed-loop control with respect toa quadratic 
performance index. In any event the optimization anal- 
ysis provides a yardstick and paves the way even if 
more classic control is desired. Further, the closed- 
loop optimization analysis indicates that all state 
variables should generally be used in the synthesis. 

Another approach to the approximately time -optimal 
closed-loop control problem is presented in Ref. 14. 
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This technique utilizes the Hamilton —Jacobi equation. 
Time-optimal reactor control is also utilized as the 
optimal control of a nuclear rocket engine if pro- 
pellant consumption is minimized.” For such a sys- 
tem, specific impulse (propellant temperature) is in- 
creased to its rated value as fast as possible while 
propellant flow is kept as low as allowable. Near the 
terminal temperature, flow is increased rapidly to 
the level for design thrust. Many constraints such as 
those due to core support thermal stresses and turbo- 
pump performance (stall and cavitation) must be con- 
sidered. Also, it is readily seen that a maximum al- 
lowable effort control is a candidate for the optimal 
reactor process as long as the cost function L is, at 
most, linear in the control. 

The practicing engineer might suggest that the 
above controls could have been synthesized by physi- 
cal arguments. Indeed, intuitive arguments were used 
to arrive at a practical suboptimal bang-bang reac- 
tivity process discussed in Ref. 15. Current optimal 
control theory puts this type of work on a firmer 
foundation and leads to a more thorough understand - 
ing of the optimal process. 

Above all, it should be realized that the engineer 
cannot drop his physical reasoning and classical con- 
trol theory to rely on mathematics alone to generate 
optimal solutions. All are necessary design ingredi- 
ents. 


QUADRATIC PERFORMANCE INDEXES 

Zero-power reactor control that minimizes aquad- 
ratic performance index involving reactivity and power 
for some desired trajectory is analyzed by dynamic 
programming in Refs. 16 and 17. The optimal control 
is found by making a nonlinear transformation of the 
neutron kinetics to a linear system discussed in Ref. 
16. Shen and Haag" applied dynamic programming to 
find a fairly complicated control process. A quadratic 
performance index involving a variation from a de- 
sired neutron level and from a desired rate of change 
of neutron level has been investigated;!’ the analysis 
utilizes the maximum principle to arrive at a fairly 
complicated optimal control. As in Ref. 16, the opti- 
mal feedback control involves time-variant gains 
(hyperbolic functions). Such solutions, although opti- 
mal with respect to a quadratic performance index, 
would be suboptimal relative to a minimal-time tra- 
jectory. 

Optimal control of a two-energy group neutronic 
model has been analyzed in Ref. 18. Here it is sug- 
gested that the optimal control really involves spatial 
variations, as well as time variations, in reactivity. 
This presents a host of challenging design problems. 


XENON POISON AND SINGULAR SOLUTIONS 

Reactor shutdown to minimize xenon peak produc - 
: : : 2 
tion has been analyzed by dynamic programming** F 
and the maximum principle.”! Power pulsing (a bang- 
bang process) was shown to produce the optimal shut- 
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down. Roberts and Smith? considered minimal-time 
reactor shutdown with xenon concentration constrained 
by some maximum allowable value. For such a pro- 
cess the optimal power variation usually includes a 
single power pulse but is complicated by a continuous 
power variation to maintain the maximum allowable 
xenon concentration if such a value is reached. Also, 
a singular solution —a state trajectory that causes the 
switching function to be zero —can be part of the op- 
timal trajectory. In other words the ideal optimal 
power usually would be zero until the xenon constraint 
value is reached, and then power is varied to maintain 
constant maximum allowable xenon until a certain 
minimal-time full-power trajectory is intercepted. 
The latter trajectory is optimal until the zero-power 
trajectory tangent to the maximum xenon line (in 
xenon —iodine state space) is reached, This new por- 
tion is then optimal to the terminal end (no poison). If 
the initial poison buildup is small, the required shut- 
down is merely zero power to the end. Both cases as- 
sume the singular condition is never reached. If the 
singular condition arises, it replaces a portion of the 
above trajectories. 

Singular solutions were analyzed in Ref. 8 but were 
not time optimal for such problems. In most studies, 
however, singular solutions have been neglected. Such 
neglect can lead to erroneous answers, as in the case 
of xenon-poison consideration for shutdown. 


Nonlinear Analysis 


Neutron kinetics ordinarily are defined by a bilinear 
system of ordinary differential equations. In other 
words the system includes a p(/)n(/) term andis linear 
in state and input (reactivity p) but not jointly linear 
in both. If reactivity is a time function, the system is 
linear and time variant. Generally, however, reactiv- 
ity is a function of neutron level (via temperature co- 
efficients and control), and the system is nonlinear. 

Although heat-exchange processes are nonlinear, 
distributed, and often time variant, it is almost al- 
ways assumed in the papers discussed below that the 
feedback systems are lumped, linear, and time invari- 
ant. Such assumptions are not always valid, but still 
the results can be usefulifone realizes the conditions. 

In the following papers, classical topological tech- 
niques’? are sometimes utilized, but most of the re- 
cent work leans heavily on Liapunov’s second (or di- 
rect) method."4 Although the second method is most 
powerful, it yields only sufficient conditions that fre- 
quently give conservative estimates of stability. The 
technique involves generation of a particular type of 
function called a Liapunov function. Fortunately some 
techniques have been devised to bypass the trial and 
error involved in the generation of Liapunov functions. 

Normally it is required that the reactor be stable — 
at least to the desired terminal state. Over certain 
portions of optimal trajectories, however, divergent 
modes might be acceptable. Such is the case for 
minimal-time startups. Only the most recent work 
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with which I am acquainted is presented here. But 
first it should be noted that significant spadework 
appeared earlier by Welton, Weinberg, Ergen, Smets, 
Gyftopoulos, Popov, Akcasu, and DeVooght, >? to 
mention a few. 

Two recent papers show that reactor stability 
without delayed neutrons does not always imply sta- 
bility with delayed neutrons. Although delayed neu- 
trons usually have a stabilizing effect, reactor engi- 
neers Should be aware of a destabilizing possibility. 
Furthermore, Smets™ shows that, if reactor -feedback 
elements exhibit a phase-lag frequency response, de- 
layed neutrons are Stabilizing. 

Because the second method of Liapunov yields suffi- 
cient but not necessary conditions, the conclusions 
are often so conservative as to have little practical 
value. This is the case despite the fact that research- 
ers continuously are reporting the latest, most non- 
restrictive criteria such as have been claimed justly 
in Refs. 35 and 36. Also, Schultz®’ has shown that 
Luré’s method yields smaller stability regions for 
reactors (and is thus less useful) than does the 
variable-gradient method. The second method does 
allow a start for detailed computer studies in the se- 
lection of certain control-system parameters. 

A useful detailed comparison of the Luré —Welton 
criterion, *!-33 Popov’s criterion,*® Aizerman’s meth- 
od,**-#° Rosen’s method,*! a topological scheme, and a 
direct application of Liapunov’s second method is 
presented in Ref. 42. Although a combination of tech- 
niques should usually be considered in practice, for 
the typical case shown in Ref. 42, a Liapunov function 
can be found directly that ensures stability over a 
Significantly larger region of operation than do the 
other methods. 

In some cases, analyses have been too specific; 
erroneous conclusions have resulted. For example, the 
conclusions reached in Ref. 34 indicate that stability 
for small perturbations without delayed neutrons does 
not imply stability for small variations with delayed 
neutrons as previously indicated.* Also, it has been 
shown in Ref. 30 that results reached in Ref. 26 are 
erroneous. That is, a reactor can even be asymptot- 
ically stable without delayed neutrons and unstable 
with delayed neutrons. 

Another interesting conclusion is formulated by 
Lellouche**: for systems with more than one reactiv- 
ity coefficient and with a prompt negative coefficient, 
a delayed positive coefficient is more stabilizing than 
a delayed negative coefficient. Also, Akcasu and 
Noble“! show that, if a prompt negative reactivity ex- 
ists, the reactor is stable if total reactivity-feedback 
transfer function (with respect to neutron level) has 
no positive real zeros and if the time integration of 
the nonprompt feedback-impulse response is negative 
for all positive time. 

Chernick** and Shotkin*® have made worthy contribu - 
tions to nonlinear reactor theory for some time. In 
fact, Chernick*’ was probably one of the first to indi- 
cate that delayed neutrons do not always provide 
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damping. In a more recent paper, Garabedian and 
Lynch’® analyze a two-energy-group multiregion 
neutron-kinetics model. Also, stability of coupled- 
core reactors (which might be economically feasible 
for nuclear rocket propulsion) with delay times is 
analyzed by Liapunov’s second method in Ref. 49. 
Such coupled systems seem to want to work together. 


Conclusions 


Classical assumptions that are usually made (some- 
times erroneously) include one or more of the follow- 
ing: zero-power reactor model, prompt-neutron kinet- 
ics, single-precursor model, linear time-invariant 
heat exchange, infinitely fast actuators, constant in- 
puts, and point (space) dynamics. In the optimal con- 
trol papers, singular solutions have been neglected in 
most cases. Also, optimality is usually assumed when 
only extremal solutions have been derived. These as- 
sumptions are sometimes valid, but in most cases 
many of them are invalid. 

Undoubtedly, trial-and-error techniques with digi- 
tal- and analog-computer simulations and classical 
linear-systems theory will continue to be useful with 
increasing help from modern techniques. Although we 
cannot hope to find analytical solutions to most prac- 
tical reactor-control problems, intuitive arguments 
would often bridge the gap between theory and prac- 
tice. This has been especially useful in showing that 
extremal solutions are sometimes indeed optimal.® 

Because of the urgency involved, many studies of 
modern reactor-control theory have involved aero- 
space or military applications. On the other hand, it 
should be realized that the same sort of analysis is 
feasible for commercial reactor systems. For exam- 
ple, minimal-time control and xenon-poison consider - 
ations are significant for commercial reactors if 
maximum operational time for a given cost is of in- 
terest. Similarly, it is important to understand the ef - 
fect of certain design parameters on stability of the 
nonlinear reactor, for safety reasons if nothing else. 
Undoubtedly there is a need for studies of cost effec - 
tiveness and optimization on commercial reactors and 
nonlinear dynamical analyses of more realistic com- 
mercial reactor models. 


* ” * 


The author is grateful to the Sandia Corporation for 
support of research in this area. Also, some of the 
ideas presented here are based on the author's pre- 
vious work at the Los Alamos Scientific Laboratory. 
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FLUID AND THERMAL TECHNOLOGY 


Improved Design Predictions 
Through Liquid-Metal Heat-Transfer Research 


By Ralph P. Stein 


One of the more frequently occurring objectives 
of basic research in heat transfer istofind prediction 
methods that will be useful for the design and analysis 
of systems in which there is technological interest. 
For convection heat transfer, the pertinent funda- 
mentals are fairly well developed. These fundamen- 
tals include the familiar conservation laws related 
to mass, energy, and momentum together with re- 
lations describing their local rates of transport. 
Direct application of these fundamentals to specific 
situations usually results in complicated mathemati- 
cal formulations that, for engineering purposes, would 
better be avoided if possible. As a result, research 
in this area has centered about application of these 
fundamentals to general situations so as to discover 
prediction methods that will be more convenient for 
engineering use. 

To use a deceptively simple and certainly familiar 
example, consider the relation 

q = hit, —t) (1) 

where q is the local rate of heat flow per unit area 
of and from a surface at temperature /, to a fluid 
with a suitably defined average temperature ¢. This 
equation defines the heat-transfer coefficient i. When 
h can be considered a constant, or treated as some 
appropriate average value, use of this equation is 
much simpler than the use of the mathematical 
formulation of fundamentals it replaces. Basic re- 
search studies discovered best definitions of ¢ for 
several general classes of applications. The heat- 
transfer coefficient could then be considered inde- 
pendent of q and (¢;—¢) and easily related to other 
variables without destroying the simplicity of Eq. 1, 
thereby supporting its practical usefulness. 

Similar potentially useful results of basie re- 
search pertaining to liquid-metal heat transfer will 
be reviewed here. These include improved methods 
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for predicting surface temperatures in reactor cool- 
ant channels and heat-transfer rates in heat ex- 
changers. They derive from applications of funda- 
mentals to analytical investigations of liquid-metal 
forced-convection heat transfer.'” 


COEFFICIENTS FOR LIQUID METALS 


With turbulent forced-convection heat transfer of 
nonmetallic fluids in ducts, the “heat-transfer-co- 
efficient concept,’ —namely, that it is possible and 
beneficial to consider gq proportional to (é, —¢) — 
has proved to be both accurate and convenient in 
most situations of engineering interest. With liquid 
metals, however, the situation is more complicated 
because, unlike the situation with nonmetallic fluids, 
heat-transfer coefficients are affected significantly 
by the actual heat-transfer boundary conditions, such 
as heat-flux or temperature distributions at the duct 
surfaces. 

Although the sensitivity of liquid-metal heat-trans- 
fer coefficients to boundary conditions is well known, 
satisfactory prediction methods that account for this 
sensitivity with engineering computations have not 
been firmly established. The implication of current 
practice is that Eq. 1 is still to be used, but with the 
complication that the heat-transfer coefficients must 
be appropriate to the specific boundary condition of 
the application. However, the many possible boundary 
conditions of interest cause much of the convenience 
and apparent simplicity of Eq. 1 to be lost. Further, 
in many cases determination of appropriate heat- 
transfer coefficients—for example, by application 
of fundamentals or experiment—is equivalent to 
“solving” the problem, after which application of 
Eq. 1 is no longer necessary. Thus basic research in 
liquid-metal convection heat transfer would search 
for prediction methods that can account for the effect 
of boundary conditions, as well as for other possible 
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effects, in more convenient forms than implied by a 
consistent use of Eq. 1. 


HOW IMPORTANT ? 


How important is this type of basic research to 
current liquid-metal technology? In many design 
applications of Eq. 1 to systems involving liquid- 
metal heat transfer, the effects of boundary conditions 
on appropriate values of / are mostly ignored. This 
can be justified by expediency —why bother compli- 
cating a situation in which, when the entire system 
is considered, other aspects may be more important 
and cannot be accounted for accurately ? 

For example, because of the excellent thermal 
properties of liquid metals, predictions of tempera- 
tures and heat-transfer rates for most design prob- 
lems in existing or proposed liquid-metal-cooled 
reactor systems are not greatly influenced by the 
accuracy with which liquid-metal heat transfer can 
be predicted. In typical liquid-metal-cooled reactors, 
because of the relatively poor thermal properties of 
present fuel materials, temperature differences be- 
tween the surface of a coolant channel and the center 
of the fuel material are roughly an order of magni- 
tude larger than temperature differences associated 
with convection of the coolant. Moreover, the perti- 
nent physical properties of these fuel materials are 
not known very accurately. In current liquid-metal 
heat exchangers and liquid-metal-heated steam gen- 
erators, structural and safety considerations usually 
result in large thermal resistance to heat transfer 
through the exchanger walls. As a result, inaccurate 
values of fh for application of Eq. 1 have little effect 
on the overall design. Also, with steam generators, 
steam-side heat-transfer predictions are highly un- 
certain. 

However, with the interest in liquid metals as in- 
dustrial coolants? as well as working fluids in fast 
reactors and in space power plants, refinements in 
design soon may require improved accuracy in pre- 
dicting liquid-metal heat transfer. How soon will 
depend on the future direction and interest of the 
related technology. 


Improved Prediction Methods 


Three examples of the results of basic research 
in liquid-metal heat transfer illustrate potentially 
useful design prediction methods. The first example 
is related to appropriate forms of Eq. 1 when applied 
to fuel-coolant channels idealized as annular spaces 
with unequal heat fluxes on opposite sides; it has its 
origins in a paper by Seban‘ in 1950. It was also 
discussed in the 1952 edition of the Liquid Metals 
Handbook® and reformulated in somewhat more gen- 
eralized terms® in 1956. Even so, there are indica- 
tions that communication of the results of basic re- 
search have not been effectively communicated for 
this particular example. This is evidenced by the 


contents of several recent surveys of liquid-metal 
heat transfer for design. 

The second example suggests an extension of 
Eq. 1 that accounts for the effect of nonuniform axial 
heat-flux distributions. The third example refers to 
improvements in predictions of heat-transfer rates 
in liquid-metal-to-liquid-metal heat exchangers for 
which indiscriminate applications of Eq. 1 can lead 
to inaccurate designs. 


UNIFORMLY HEATED ANNULAR SPACES 


Let subscripts 1 and 2 identify quantities associated 
with the inner and outer walls of an annular space. 
Then Eq. 1 as applied to the inner wall can be written 
as 


Gy = My (ts 4 — 2) (2) 


with ¢ the bulk or heat-content mean temperature of 
the fluid. When q; and q2 are known, the bulk tem- 
perature is computed easily by simple heat balance. 
But what is the proper value for the heat-transfer 
coefficient ;, so that the surface temperature /. ; can 
be predicted? For liquid metals, /, is influenced® 
strongly by the ratio of the heat fluxes q,/q,. Let 
this ratio be represented by y, and consider the heat- 
transfer coefficient to be a function of y —that is, 
hy(y). Then Ref. 6 shows that 


h,(0)/hy(vy) =1l=a,y (3) 


where a; is a positive number which depends on the 
Peclet and Prandtl numbers of the fluid, onthe diame- 
ter ratio of the annular space, and, when the heat 
transfer is not “fully developed,” on duct axial posi- 
tion also but which is independent of y. Using Eq. 3, 
we can rewrite Eq. 2 as 


G1 — 44q2 = Ny (0) (ts 4 — t) (4) 


from which, with a; and /2,(0) known, leit can be com- 
puted easily. For fully developed heat transfer, /2;(0), 
the heat-transfer coefficient for an internally heated 
annulus can be obtained from existing correlations’; 
a, can be obtained from the detailed computations of 
Kays and Leung® over wide ranges of Peclet and 
Prandtl numbers and diameter ratios.' Alternately, 
a, can be obtained from existing correlations for 
fully developed heat-transfer coefficients in sym- 
metrically heated annular spaces—that is, from 
h,(1)—since 


ay = [hy(1) — 2, (0)] /h,(1) (5) 


from Eq. 3. By use of Eq. 5, another form of Eq. 4 
can be written as 


(44 —42)/24(O) + ¢2/hy(1) = 4,4 — 4 (6) 
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which, with existing correlations for fully developed 
values of /2,(0) and (1), may be the most convenient 
form for predicting /, , for design purposes. 

The appropriate form of Eq. 3 for the outer wall 
of the annular space is 


hy(0)/hg(y) = 1 —agy (7) 


where now \=4;/q,. Similar algebraic manipulation 
leads to 


(42 — 1 4)/N2(0) + G4/h2(1) = b2 — 4 (8) 


for predicting the surface temperature of the outer 
wall of the annular space, /, ». Correlations of fully 
developed values for /t,(0) and /,(1) are also avail- 
able.':* Note that Eqs. 6 and 8 give simple expres- 
sions for predicting /, 4 when q; = 0 and for /; » when 
qd = 9. 

Relations similar to these are given in Refs. 5 
and 10. 


NONUNIFORM AXIAL HEAT FLUX 


Equation 1—-as applied to liquid-metal heat trans- 
fer with q a function of axial position / measured 
from the duct inlet—results in/ being a function of 
the flux distribution and of the axial position (except 
for certain special cases). In fact, for simple “un- 
chopped” cosine flux distributions, /; —when defined 
by Eq. 1—will always be zero at the fuel-channel 
outlet. As a result the customary definition of h as 
given by Eq. 1, with / the usual bulk temperature, 
becomes cumbersome to use because then /: must be 
determined in advance as a function of the flux distri- 
bution and of the axial position. 

A recently suggested’ extension of Eq. 1 avoids 
this difficulty but involves the use of higher ordered 
dimensionless thermal-resistance coefficients R,, 
associated with successive derivatives of q with re- 
spect to 7. The coefficients R, are functions of 
Peclet (or Reynolds) and Prandtl numbers and duct 
shape but are independent of the flux distribution 
and, for sufficiently large /, are also independent of 
axial position. This extension of Eq. 1 for sufficiently 
large / can be written as 


qt 2 Ry (i"q/dx") = h(t, —t) (9) 


where x=/1/D (D is the duct hydraulic equivalent 
diameter) and h is the fully developed heat-transfer 
coefficient for uniform heat-flux distribution for a 
circular tube, an internally or externally heated 
annular space, or an annular space with equal heat 
fluxes from both sides. Equation 9 is a dimensional 
form of Eq. 49 in Ref. 1. Note that when q is uni- 
form —that is, independent of x—-Eq. 9 reduces to 
Hq. 1. 
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The coefficients R, decrease rapidly with increas- 
ing n, so that in most cases only a few terms of the 
infinite series need be retained. In other cases the 
higher order derivatives may be zero. For example, a 
reasonable approximation of q(x) for many reactor 
applications is 

q = (4d 4.) (D/L)(L/D-%), (10) 
where L is the fuel-coolant-channel length. Then 
Eq. 9 is simple to use to predict maximum fuel or 
channel surface temperatures because only the co- 
efficients R, and R, are required. 

For the coolant channel idealized as the space 
between infinitely wide parallel planes, the coeffi- 
cients R,, can be computed from 


R, = (-1)" (Nu/2) (Pe/4)" 5 C,/x2" (11) 
k=1 


where Nu is the appropriate fully developed Nusselt 
number corresponding to the fully developed heat- 
transfer coefficient for uniform heat-flux distribution 
h used in Eq. 9 and Pe is the Peclet number. Repre- 
sentative values of C, and A} are given in Tables 1 
and 2 as functions of Reynolds and Prandtl numbers. 
These values were computed from the “Fundamental 
Eigenvalues Problem”! related to liquid-metal heat 
transfer in heated parallel-plane ducts.!! The com- 
putations were based on use of von Karman’s uni- 
versal velocity profile with eddy diffusivities for 


Table 1 VALUES OF C, AS A FUNCTION OF REYNOLDS 
AND PRANDTL NUMBERS 











Prandtl Reynolds Cy for Bq, 14 

number number C, C, Cy Cy Cc; 
0.003 30,000 0.1798 0.0509 0.0232 0.0134 0.0086 
100,000 0.1619 0.0485 0.0219 0.0129 0.0082 
300,000 0.1280 0.0413 0.0183 0.0110 0.0070 
1,000,000 0.0802 0.0269 0.0119 0.0075 0.0049 
0.01 30,000 0.1538 0.0463 0.0205 0.0117 0.0072 
100,000 0.1181 0.0388 0.0170 0.0102 0.0065 
300,000 0.0753 0.0255 0.0113 0.0072 0.0047 
1,000,000 0.0340 0.0199 0.0056 0.0038 0.0026 
0.03 30,000 0.1127 0.0363 0.0149 0.0082 0.0045 
100,000 0.0694 0.0236 0.0103 0.0064 0.0040 
300,000 0.0336 0.0119 0.0055 0.0037 0.0025 
1,000,000 0.0120 0.0046 0.0023 0.0017 0.0012 





Table 2 VALUES OF Aj AS A FUNCTION OF REYNOLDS 
AND PRANDTL NUMBERS 





2 
Prandtl Reynolds Ay for Eq. 11 eae 
number number yi Ni MR N N 








0,003 30,000 11.67 44,92 100.84 178.42 278.56 
100,000 13.16 51.24 116.08 206.54 322.96 

300,000 16.91 68.66 158.91 284.79 446.08 

1,000,000 26.63 121.16 290.52 523.80 819.36 


0.01 30,000 13.87 53.84 121.96 216.96 339.36 
100,000 18.42 75.34 175.05 314,10 491.98 

300,000 28.18 129.58 311.83 562.13 879.13 

1,000,000 52.75 295.50 733.46 1320.98 2056.17 


0.03 30,000 19.20 77.79 180.43 323.34 506.69 
100,000 30.38 140.43 338.16 609.74 953.43 

300,000 53.32 297.39 736.65 1327.76 20€7,21 

1,000,000 116.24 785.12 1977.85 3551.05 5501.29 
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heat transfer taken equal to those for momentum 
transfer and assumed to be uniform and equal to 
their maximum values at the duct center. Computa- 
tions by other methods will not give significantly 
different results, except possibly for large Peclet 
numbers where very accurate values of the co- 
efficients R, are not needed. 

Reference 1 shows that, for cosine heat-flux distri- 
butions, when 


(7/4) (PeD/L) « x% (12) 

Eq. 9 reduces to the very simple form 
q + Ry dq/dx) = h(t, —t) (13) 

Further computations” indicated that 
(1/4) (PeD/L) = 0.05X% (14) 


is usually sufficient to justify the use of Eq. 13. 
Applications of Eq. 13 to predict maximum fuel and 
surface temperatures for typical liquid-metal-cooled 
fast reactor designs revealed that neglecting the 
second term on the left side of Eq. 13 gave values 
that were less than 10°F too low. Thus, for practical 
design purposes, and with / taken equal to fully 
developed heat-transfer coefficients for uniform flux 
distributions, Eq. 13 will be sufficiently accurate for 
application to cosine heat-flux distributions when 
Eq. 14 is satisfied. Accurate predictions at the 
channel outlet, or in general with larger heat fluxes 
than those allowed in current designs, will require 
full use of Eq. 18. 

Reference 1 also discusses the cosine heat-flux 
case when Eq. 14 is not satisfied, generalizations of 
Eqs. 6 and 8 for nonuniform heat-flux distributions, 
and cases for which the heat transfer is not fully 
developed. Tables 1 and 2 give the values of C, and 
\? (which were not available when Ref. 1 was written); 
they will be found useful for analyses of these special 
cases as applied to fuel-coolant channels idealized 
as the space between infinitely wide parallel planes. 


PREDIC TING RATES IN HEAT EXCHANGERS 


The necessity of obtaining alternate or extended 
forms of Eq. 1 in the previous two examples is re- 
lated to the sensitivity of liquid-metal heat-transfer 
coefficients —as defined by Eq. 1—to the duct heat- 
flux boundary conditions. With heat exchangers, these 
boundary conditions depend on the operating condi- 
tions of the exchanger and, except for special limiting 
cases, cannot be specified in advance. The situation 
is further complicated with liquid metals because a 
significant portion of the overall heat-transfer rate 
occurs along the entrance regions of the exchanger 
where heat-transfer coefficients are very large and 
vary with duct axial position—i.e., along regions 
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where the heat transfer is not fully developed. These 
considerations suggest that the design of liquid- 
metal-to-liquid-metal heat exchangers based on the 
traditional applications of Eq. 1, as used for the 
turbulent flow of nonmetallic fluids, can lead to 
important inaccuracies. 

Current basic research on simple double-pipe heat 
exchangers!” has shown that traditional methods of 
heat-exchanger design, when applied to heat transfer 
between liquid metals, leads to inaccurate predictions 
in certain operating ranges —although in other ranges 
inaccuracies are not too important. When these in- 
accuracies are serious, oversized heat exchangers 
result. 

The practical significance of this research is seen 
in Tables 3 and 4, which compare predictions of 
concentric-tube countercurrent-flow heat-exchanger 
lengths (or heat-transfer areas) required for various 
specified overall heat-transfer rates. Heat-exchanger 


Table 3 ERRORS IN PREDICTING 
HEAT-EXCHANGER* LENGTHS 


Error in predicted length, 





Specified 

efficiency H=0.1 H=05 H=1.0 
0.3 202 56.9 
0.4 178 39.2 
0.5 167 65.0 27.3 
0.6 154 55.4 16.9 
0.7 143 46.0 11.4 
0.8 136 36.9 8.0 


0.9 122 30.9 2.6 


*Countercurrent flow; K = 0.1, A,,= 0.05. 


Table 4 ERRORS IN PREDICTING 
HEAT-EXCHANGER* LENGTHS 


Errors in predicted length, % 








Specified 

efficiency K, =0 K,, = 9.05 K,, = 0.1 
0.3 206 
0.4 155 
0.5 115 65 46 
0.€ 87 55.4 38.1 
0.7 71.2 46.0 31,2 
0.8 5.6 36.9 25.8 
0.9 03.2 30.9 0.0 
*Countercurrent flow; K = 0.1, H = 0.5. 


lengths were computed by simple applications of 
Eq. 1, with 2” assumed to be independent of axial 
position and equal to values corresponding to the 
boundary condition of uniform heat flux; these are 
compared with accurate predictions by a detailed 
analysis based on fundamentals.!? The detailed analy - 
ses show that the operating conditions of the ex- 
changer are characterized by the threedimensionless 
parameters H, K, and K,, so that these parameters 
must be defined before the comparisons can be 
understood. 

The parameter H is the heat-capacity flow-rate 
ratio defined by 


H = c,W,/c,W, (15) 
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where c is the fluid specific heat, W is the mass 
rate of flow, and subscripts 1 and 2 refer to the 
tube and annulus sides of the exchanger. The parame- 
ter K is the relative thermal resistance of the fluid 
in the annular space and is defined as 


K = (kykt/kk})(1 —R)/R (16) 


where k is the fluid thermal conductivity, R is the 
diameter ratio of the annular space, and k* accounts 
for the contribution of turbulent flow to the heat 
transfer with satisfactory accuracy when Peclet 
numbers are less than 1000.' For the tube side of 
the exchanger 


k, ~ 0.875 + 0.00312(Pe,)":® (17) 


For the annulus side with R close to unity —that is, 
when the annular space is narrow, 


kt} = 0.965 + 0.00333(Pe,)°-® (18) 


The parameter XK, is the relative thermal resistance 
of the tube wall which, when the wall is thin, is de- 
fined by 

Ky = (kyRT/R)(b/7,) 0/7, «1 (19) 
where k, is the thermal conductivity of the wall 
material, ) is the wall thickness, and 7; is the tube 
radius. The computations assume R is close to unity, 
and thus K will be a very smallnumber. The parame- 
ter K, will usually be a small number also. A value 
of K of 0.1 and a range of values of K, from 0 to 0.1 
were chosen as reasonable representations of liquid- 
metal-to-liquid-metal concentric-tube exchangers 
with narrow annular spaces. 

In Tables 3 and 4 the specified overall heat- 
transfer rate is shown as an “efficiency,” defined as 
the ratio of the actual heat-transfer rate to the rate 
obtainable for an exchanger with infinite heat-trans- 
fer area at the same operating conditions. The length 
(or heat-transfer area) comparisons, shown as per- 
cent errors of the length predicted by use of Eq. 1 for 
three values of H, are always larger than the lengths 
computed by the accurate detailed analysis. The 
errors are extremely large for values of H of 0.1 but 
negligibly small for H=1 with moderate to high 
efficiencies. 

Errors in predicted lengths result from two effects: 
(1) the sensitivity of liquid-metal heat-transfer co- 
efficients to heat-flux boundary conditions and (2) the 
high rates of heat transfer in the entrance regions 
of the exchanger. For the not unusual case of H = 1 
and a sufficiently long exchanger (equivalent to a high 
efficiency), heat fluxes along the central regions of 
the exchanger become independent of axial position — 
so fully developed heat-transfer coefficients cor- 
respond to the boundary condition of uniform heat 
flux, as used when applying Eq. 1; thus errors are 
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small. The increase in percent errors of predicted 
length with decreasing efficiency for the case of 
H = 1 reflects the importance of the entrance regions 
of the exchanger where heat fluxes are significantly 
higher than those in the central regions. For values 
of K, significantly larger than 0.1, the wall thermal 
resistance begins to control the local heat-transfer 
rates; as a result the local heat-flux distributions 
and overall heat-transfer coefficient approach uni- 
formity, so that the traditional applications of Eq. 1 
become sufficiently accurate. 

The computations of the detailed analysis based on 
fundamentals are complex and require the use of a 
large digital computer.” Simpler computational pro- 
cedures are preferred for engineering applications, 
especially when heat-exchanger performance must be 
predicted over large ranges of conditions when the 
operation of a complete power-plant cycle is ex- 
plored. Improved design and performance relations, 
based on extensions of the familiar effectiveness (or 
efficiency) vs. NTU types of equations,'® have been 
suggested as part of the basic research investiga- 
tions.'’? These extensions introduce two correction 
factors: n, which is the ratio of the actual fully de- 
veloped overall heat-transfer coefficient to the fully 
developed overall heat-transfer coefficient based on 
uniform heat-flux boundary conditions, and ¢, called 
the “effectiveness coefficient,” which accounts for the 
high local heat-transfer rates in the entrance regions 
of the exchanger by multiplying the exponential terms 
in the usual effectiveness vs. NTU equations. 

Both n and @ depend on H, K, and kK, , but @ also 
depends on the heat-exchanger length (or, equiva- 
lently, the heat-exchanger efficiency). However, for 
sufficiently large efficiencies, ¢ becomes independent 
of length and is called the “fully developed effective- 
ness coefficient.” Table 5, an extension of a similar 
table in Ref. 2, shows typical values of these correc- 
tion factors for the circular-tube countercurrent flow 
heat exchanger with a narrow annular space. The 
“fully developed effectiveness coefficient” is denoted 
by $(%) and applies for values of the heat-exchanger 


Table 5 CORRECTION FACTORS FOR THE CONCENTRIC-TUBE 
COUNTERCURRENT-FLOW HEAT EXCHANGER 





K,=09 K, = 9.05 Ky = 9.1 
K H o(~) n o(~) 7 o(~) uf 
0.01 0.1 0.6 4.74 0.8 2.39 0.8 1,851 
0.5 0.78 1.468 0.85 1.229 0.88 1,160 
0.9 0.966 1.037 0.976 1.025 0.982 1.018 
1.1 0.971 0.971 0.979 0.979 0.985 0,985 
2.0 0.85 0.851 0,881 0.890 0,91 0.916 
0.05 0.1 0.63 3.50 0.8 2.18 0.82 1,767 
0.5 0.792 1,318 0.856 1.206 0.892 1,147 
0.9 0.968 1.033 0.978 1.023 0.983 1.017 
1 0.973 0.974 0.981 0.981 0.986 0,986 
2.0 0.86 0.863 0.888 0.898 0.914 0.922 
10.0 0.7 0.770 0.79 0.823 0.83 0.861 
0,10 0.1 0.67 2.76 0.8 1.990 0.827 1.678 
0.5 0.811 1.269 0.867 1.181 0.899 1.133 
0.9 0.971 1,029 0.979 1.021 0,984 1.016 
1 | 0.975 0.977 0.982 0.983 0.986 0.987 
2.0 0.86 0.878 0,894 0.908 0.918 0.929 


10.0 0.75 0.795 0.79 0.842 0.83 0.876 
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efficiency larger than 0.7. The values of $(%) may 
also apply for much smaller efficiencies, but further 
investigations are required to verify this possibility.’ 
Most of the results shown in Tables 3 and 4 can be 
reproduced with satisfactory accuracy by the values 
of » and $(*) in Table 5 and by noting that 7 = 1 and 
@(%) = 1 correspond to the traditional applications of 
Eq. 1 for heat-exchanger design as used with non- 
metallic fluids in turbulent flow. 


CONCLUDING REMARKS 


These three examples of contributions of basic 
research to design are based on fundamental studies 
of relatively simple models. However, actual liquid- 
metal-cooled reactor fuel channels are not usually 
simple annular spaces or circular tubes —and actual 
heat exchangers rarely have simple double-pipe con- 
figurations. Even so, the research suggests general 
principles that often can be made applicable to more 
complex systems. For example, procedures for ap- 
plying the results of basic research on simple annular 
spaces to fuel-element rod bundles or heat-exchanger 
tube-in-shell configurations have been suggested"! 


and have been,'® and continue to be, tested by ex- 
periment. 
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OPERATING EXPERIENCE 


NPD Performance Indicates Merit 
of On-Power Refueling 


By Karl K. Brown 


With enriched-uranium fuel, reactors can be de- 
signed to be refueled during plant shutdowns that 
coincide with annual turbine-generator overhauls. 
But when the fuel is natural uranium, power reac- 
tors must be refueled at least every 6 months to 
maintain reactivity. In Canada, where D,O-moder- 
ated natural-uranium-fueled pressure-tube reactors! 
are favored, there is strong economic motivation to 
refuel while the reactor is operating at power so as 
to achieve high burnup and low fuel cost (Table 1). 
With successful on-line refueling comes the bonus 
of higher plant availability. 

Canadian experience in refueling the 20-Mw(e) 
NPD plant (abridged in Ref. 2) is significant not 
only to their plants and to competing systems but 
also to related designs such as the U. S. Heavy Water 
Organic Cooled Reactor’ (HWOCR) which also are to 


mh 1 








have on-line refueling. All committed Canadian nu- 
clear power plants —the two 500-Mw(e) CANDU type 
units at Pickering, Ontario, scheduled for operation 
in 1970 and 1971, and the 250-Mw(e) CANDU-BLW, 
also scheduled for 1971 —vwill refuel at power. 
On-power refueling machines are already installed 
in the new 200-Mw(e) CANDU reactor at Douglas 
Point, in the 20-Mw(e) NPD (Figs. 1 and 2), and also 
in a research reactor—the NRU. Operation of the 
NRU on-power refueling equipment began in 1957, 
although serious design consideration dates back to 
1949. The more than 2 years of experience with on- 
power refueling machines on the NPD power reac- 
tor are statistically inadequate, but the prognosis 
is for high unit availability and low forced-outage 
rates. Plant unavailability resulting from refueling- 
machine problems at NPD amounted to only 2.44% 





Fig. 1 NPD on-power refueling machine during commissioning in 1962. 
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Why Refuel at Power? 


Natural uranium is an essential feature in the Canadian approach to 
nuclear power. In effect, they choose to enrich the moderator rather 
than the fuel.4 With this choice, howevei, reactivity requirements com- 
mit them to refuel the reactor no less frequently than semiannually. 
Utility energy demands, the cost of shutdowns, and fueling costs make 
it desirable to refuel while the reactor continues producing power. 

Table 1 shows comparative power costs for the first pair of 500- 
Mw(e) units in the Pickering generating station. These estimates are for 
conditions in Canada and are based on long-term Ontario Hydro finance 
criteria. Capital and operating costs of on-power refueling are included 
and are based on a burnup of 223 Mw-hr per kilogram of uranium, a net 
unit efficiency of 29.7%, and a unit fuel cost of (Canadian) $50 per kilo- 
gram of uranium. The low fuel cost of 0.75 mill/kw-hr can only be 
achieved with on-power refueling. In the Ontario Hydro power system, 


, the difference in incremental operating costs between coal and nuclear 
generation is ~2 mills/kw-hr in favor of nuclear with on-power re- 
fueling. 


Table 1 ESTIMATED POWER COST* 
WITH REFUELING ON-POWER VS. OF F-POWER 


With off-power refueling 














With aes 
on-power Two times Three times Four times 
Category refueling per year per year per year 

Capital charges 2.18 2.24 2.30 2.36 
Operations and 

maintenance 0.42 0.43 0.44 0.45 
Fuel 0.75 2.10 1.50 1.40 

Total 3.35 4.77 4.24 4.21 

Relative 100% 143% 127% 126% 


(Reference) 





*Expressed as Canadian mills per kilowatt-hour and assuming 5.25% 
interest, 30-year life, 80% capactiy factor, and 6.69% fixed charges for 
the Pickering generating station with two 500-Mw(e) net units. 


If off-power refueling is employed, the reactivity-limited achievable 
burnup raises fuel cost and, thus, power cost. For instance, with semi- 
annual off-power refueling, which is barely possible from a reactivity 
point of view, burnup is only 79 Mw-hr per kilogram of uranium and the 
complete core change triples the fuel cost from 0.75 mill/kw-hr to 2.1 
mills/kw-hr. With off-power refueling three times annually, burnup can 
reach 110 Mw-hr per kilogram of uranium and the fuel cost is only ap- 
proximately doubled to 1.5 mills/kw-hr. Quarterly refueling, which is 
unacceptable owing to seasonal peak-power requirements, would in- 
crease burnup to 120 Mw-hr/kg. 

It can be assumed that one off-power refueling could coincide with an 
annual shutdown; others would cause shutdowns. The cost of unavail- 
ability is interdependent with the cost of fuel and also depends on both 
design concept and hardware. For infrequent off-power refueling, each 
additional shutdown reduces the annual availability by ~2% —which 
effectively increases the energy unit cost due to capital and operation 
(fixed charges) by 2% for each shutdown. 

Cost of the entire on-power refueling system for a large station (e.g., 
Pickering) is estimated at <3% of total plant capital cost. It is unlikely 
that the cost of off-power refueling systems would be less. In fact, to 
secure high reactor availability, a faster machine would be desired; so 
machine cost might actually increase. 











during the 26-month period ending Mar. i, 1966. machine port and the machine automatically posi- 

As refueling-machine performance improves, fuel- tions it in the machine “head,” the machine tra- 
scheduling experience should also improve. verses to the reactor, clamps (leaktight) onto the 

) proper pressure-tube end fitting (Fig. 4), removes 
NPD Refueling Concept the seal plug (Fig. 5), and sets it aside within the 

head. Concurrently, a discharge machine similarly 

The NPD overall refueling system is shown in Fig. clamps onto the end fitting at the other end of the re- 


3. After an operator loads a new fuel bundle into the actor pressure tube and removes the seal plug there. 
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Fig. 2 Douglas Point on-power refueling machine during recent commissioning. This and later designs 
provide shutdown shielding between reactor and refueling machines. 
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Fig. 3 NPD refueling system. After the machines clamp onto end fittings at both ends . 
of the reactor pressure tube and remove seal plugs, the vam in one machine pushes a 
new fuel bundle into one end of the tube so that a spent bundle is pushed out the other end ——_— ——- 
and into the other machine; then the machines replace the seal plugs. The machines are i 
operated from the reactor control room. 




















Fall 1966 


Then the charging machine rams the new fuel into the 
pressure tube, pushing a spent-fuel bundle into the 
discharge machine, the machines replace the seal 
plugs in the end fittings, and the machines unclamp. 
Adjacent pressure tubes are refueled from opposite 
ends to balance power generation within the core. 








Feeder iN) 











Fuel bundle 





Seal plug : 


Fig. 4 End fitting of stainless steel on each end of Zir- 


2 


caloy-2 reactor pressure tube has seal plug that refueling 
machine removes for access lo fuel bundles. 


\— Plug body 


~~ 
tc Flexible disk seal 
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End i 
fitting Seat 


Plug 
body 
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Coolant pressure deflects 
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sealing force 


Fig. 5 Seal plug has flexible-disk seal (stainless steel) on 
which reaclor-coolant pressure acts to tighten disk against 
sealing faces. 
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Refueling-Machine Performance 


When the refueling machines were first conceived, 
several obstacles were obvious: (1) leakage at the 
end fittings could cause high D,O losses, (2) radio- 
activated equipment would complicate maintenance, 
and (3) sticking of fuel, seal plugs, and clamps could 
require plant shutdown for repairs. How well the de- 
sign has overcome these obstacles can be seen from 
consideration of refueling-machine performance. 


COMPONENTS 


Figure 4 shows the NPD end fittings, which con- 
tain D,O at 1050 psig, 530°F. Although several hun- 
dreds of tests were performed in the laboratory on 
the NPD design, a major leak occurred in the reac- 
tor in December 1962 when an elastomer seal ex- 
truded through a sealing face because of inadequate 
clamping forces. Since correction of this deficiency, 
452 routine on-power fueling operations have been 
carried out (904 clampings) without further trouble 
(as of Mar. 1, 1966). 

Seal plugs (Fig. 5) in the 264 end fittings in NPD 
must be routinely removed and replaced during re- 
fueling without serious chronic D,O loss. This re- 
quirement is met in two ways: (1) development and 
use of a long-life self-energizing seal and (2) con- 
tinued recovery of any escaped steam or liquid from 
the reactor-vault atmosphere or floor. 


LEAKAGE 


Low-leakage seal plugs are required whether fuel 
is changed under on-power or off-power conditions. 
To prove that each plug is leaktight, the machine 
automatically tests it before the machine seal is un- 
clamped. Spare seal plugs are carried in the ma- 
chine magazines and can be installed on the pressure 
tubes while the reactor is at power. There have been 
no outages or deratings due to seal-plug leakage. 

The average chronic escape rate of heavy water 
into the reactor vault for all reasons during 1964 
and 1965 was 280 g/hr. The average recovery rate 
was 262 g/hr during the same period, so the aver- 
age loss rate was 18 g/hr, or ~1 Ib/day. 


RADIATION LEVELS 


The refueling-machine heads and associated equip- 
ment in the reactor vault are subjected to high ra- 
diation levels. Although no outages or failures have 
resulted from radiation degradation, electrical wir- 
ing, hoses, and internal-machine elastomers have 
been replaced. Neoprene hoses lined with Buna-N 
rubber are now used to achieve good thermal sta- 
bility. 


MAINTENANCE 


Maintenance of the refueling-machine heads re- 
sults in personnel exposure to radiation produced by 
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direct neutron activation of surface contaminants as 
well as fromcorrosion products, fission products, and 
tritium. Although there has never been a case where 
a person has exceeded the allowable exposure of 5 
rems/year, maintenance and modification of the re- 
fueling machines during 1964 and 1965 caused a to- 
tal of 27 man-rems/year exposure. The present goal 
is <6 man-rems/year per unit (two machines). 

Because the reactor vault is inaccessible to per- 
sonnel and the machine could stick on an end fitting, 
the design permits manual emergency operations 
through openings in a rotating shield at the end of the 
reactor vault. Two machine outages have required 
the use of these openings at NPD. A special cone- 
shaped bung was inserted manually into the reactor 
pressure tube after the machine was detached. This 
has proved to be a relatively simple and straight- 
forward operation. 

In the Douglas Point and Pickering designs, this 
problem has been further reduced by placing shut- 
down shielding between the reactor and the refueling 
machine. 

Emergency provisions were also made at NPD to 
permit a refueling-machine head containing jammed 
fuel to be transferred and immersed in the spent- 
fuel canal remotely. Although fuel has jammed in the 
head on two occasions, the emergency provisions 
proved unnecessary because radiation levels were 
much lower than indicated by the conservative cal- 
culations on which the design was based. Thus these 
events were handled with relative ease and caused 
no outages. 


NPD Refueling-System Results 


The NPD has operated with on-power refueling 
routinely since Jan. 1, 1964—the data presented 
here cover experience from then until March 1966. 
Clearly, 2 years of data for one unit during a period 
of plant immaturity are insufficient for firm conclu- 
sions. Early problems are still being identified, and 
development is in progress, so performance should 
improve. 


STATION OPERATION 


During the entire year 1964 and the first half of 
1965, NPD was run to provide base load. From July 
to November 1965, numerous planned training and 
modification outages and some forced and mainte- 
nance outages occurred. During the 1965-1966 win- 
ter peak (December to February), the unit was base 
loaded with no planned outages. 

No attempt is made here to discuss the overall 
unit statistics, nor have any attempts been made to 
adjust the records for tests, modifications, etc. The 
overall values are presented solely for the purpose 
of providing perspective for the refueling-machine 
Statistics. 
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OUTAGES 


Table 2 shows the unavailability of NPD for all 
reasons, together with unavailability caused by the 
refueling system. Although the refueling system 
caused three outages during 1964 and four during 
1965, it caused no outages during January or Feb- 
ruary 1966. The average outage lasted 66 hr. Two 
limit-switch deficiencies on the refueling-machine 
heads accounted for 75% of the lost time. 


Table 2 UNAVAILABILITY HISTORY FOR NPD STATION 





Station Caused by refueling 








unavailability nas 
for all* reasons, Unavail- No. of 
Period % ability, % outages 
1964 14,14 0.63 3 
1965 23.30 4.67 4 
1966 (January 
and February) 0.06 0.00 0 
Average per 
year 17.26 2.44 3.2 





*All includes experiments, tests, modifications, maintenance, and 
all other events that caused NPD to be unable to produce electricity. 


The forced-outage rate for the winter period isa 
vital statistic in the Ontario Hydro system and is 
particularly important in planning reserve capacity. 
The 1964-1965 2-month base-load run was during 
December and January; in the winter of 1965-1966, 
the base-load peak was during a 3-month run, De- 
cember to February. During these two runs there 
were no outages caused by the refueling system. The 
perfect performance during this period is not sta- 
tistically reliable, of course, and NPD will suffer 
forced outages caused by the refueling system in the 
future. However, by careful planning they expect to 
achieve a very low forced-outage rate during the 
system seasonal peaks. 

The achievement of 2.44% refueling-machine un- 
availability might be interpreted as adequate or ac- 
ceptable performance. However, the unavailability 
data can be, and in this case are, completely mis- 
leading. Unavailability cannot substitute for forced- 
outage rate in assessing system reserves and is a 
poor yardstick for economic evaluation. Perhaps its 
greatest virtue is the ease of measurement. 


CAPABILITY 


For economic purposes, unavailability is a poor 
yardstick because it does not include derating or en- 
ergy consumed during outages and is too awkward 
for sophisticated treatment of coincident troubles. 
Accordingly, “capability” is defined as a measure of 
ability to produce net energy. Its complement, “in- 
capability,” is a measure of unit imperfection. Where 
availability uses a time base, capability uses an en- 
ergy base. 

Table 3 represents NPD capability data during the 
on-power experience period; it also shows the plant 


need 
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Table 3 NPD STATION CAPABILITY STATISTICS 





Station factor, % 


Station incapability, % 
Caused by refueling system 
For all reasons 











Operating = Direct Indirect 
Period policy Capacity* Capabilityt Outage Derating Total Outage outage Derating Total 
1964 Base load 82.6 82.9 14.93 2.14 17.07 0.8 2.5 0.2 3.5 
1965 Modification 70.3 71.3 25.00 3.70 28.70 5.4 1.4 1.5 8.3 
and training 
1966 (January 
and February) Base load 99.9 99.9 00.07 0.00 0.07 0 0.0 0.0 0.0 
Average annual 78.2 78.9 18.40 2.70 21.10 2.87 1.80 0.79 4.05 





*Capacity factor = actual net energy produced = 100 


net energy unit capable of producing x 100 


tCapability factor = ——___—_____& _— ducing 


net energy with perfect continuous operation at full capacity’ 


net energy with perfect continuous operation at full capacity’ 


incapability caused by refueling-system troubles. 
For the latter the major difference between the aver- 
age plant unavailability, 2.44%, and the average plant 
incapability, 4.05%, corresponds to energy consumed 
during shutdown, indirect outages, startup losses, 
derating during xenon buildup after outages, and de- 
rating as a result of refueling-machine trouble. 

During 1965 the incapability rose to 8.3% for two 
reasons: (1) two limit-switch problems on the re- 
fueling machine increased direct-outage incapability 
to 5.4% (development work was required to correct 
the problems) and (2) equilibrium-fuel operation in 
1965 was more representative of normal conditions 
and caused the derating incapability caused by re- 
fueling-system trouble to increase to 1.5%. Modifi- 
cations made in late 1965 and improved maintenance 
were expected to yield a major improvement in 
refueling-system performance during 1966. 

Table 4 gives present performance targets. Al- 
though these targets are not vital to the economics 


marized here, covers the time from reactor startup 
to near-equilibrium burnup. (Equilibrium occurs 
when core-average burnup reaches a constant value 
with specified burnup profiles.) Figure 6, the actual 
core-average burnup, shows that NPD is approaching 
equilibrium — although the burnup profiles are not 
yet quite as had been specified. At equilibrium, twice 
the asymptotic value is the discharge burnup; thus 
the extrapolated value from this curve will be ~10% 
over the design discharge burnup (152 Mw-hr per 
kilogram of uranium). 

In the ideal natural-fuel schedule, fuel is simply 
shifted from channel to channel until equilibrium 
(1200 x 10° Mw-hr) is reached. Then new fuel is 
added at the equilibrium rate (Fig. 7), and there is 
no further shifting of bundles among channels. How- 
ever, the ideal schedule cannot be achieved in prac- 
tice because of physical limitations of the system — 
mainly fuel movement, maximum burnup, and flux 
distribution. Furthermore, NPD scheduling was also 


Table 4 PERFORMANCE GOALS FOR MATURE STATIONS AND REFUELING SYSTEMS 





Station performance for all reasons 


Lifetime Winter-peak 


Station performance caused by refueling system 


Winter-peak 











capacity forced-outage mene. 5 forced-outage Oe 

Station factor, % rate, % Outage Derating Total rate, % Outage Derating Total 
NPD-2 1.50 0.33 0.85 0.15 1,00 
Douglas Point-1 80 3.30 6.50 7.20 0.56 1,11 0.18 1,29 
Pickering-1 80 4.70 9.00 10.00 0.55 1,06 0.25 1,31 





of this concept, they are desirable to minimize en- 
ergy cost. The real balance is reached where in- 
creased performance is offset by increased capital 
and operating costs for each variable. McConnell 
concludes’ that “although we are very pleased with 
the average fueling-system incapability experience 
of 4.05%, we are confident of achieving 1% in NPD.” 


Fuel-Scheduling Experience at NPD 


The payoff for on-power refueling with the equip- 
ment described comes when it performs to help 
achieve design burnups on an equilibrium basis and 
to increase station availability. Fuel-scheduling ex- 
perience so far, as reported by Fanjoy® and sum- 


complicated by our use of depleted fuel in the core 
center to compensate for the excess reactivity caused 
by the new fuel in the initial loading. Thus Fig. 7 also 
shows the schedule calculated® for adding new fuel 
bundles on the basis of the depleted-fuel effect and 
practical shift rules. (Fuel bundles are kept within 
the design burnup by shifting the bundles removed 
with low burnups to alternate channels and discharg- 
ing the others.) 


THE ACTUAL SCHEDULE 


Figure 7 also shows the actual numbers of bundles 
shifted between channels and the number replaced by 
new bundles. Decisions on whether and where to move 
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Fig. 6 Actual core-average burnup is approaching equi- 
librium value; discharge burnup is twice the asymptotic 
value. 
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Fig.7 Refueling schedules—theoretical ideal, calculated ® 
and actual—during the time from reactor startup with all- 
new fuel to the achievement of near-equilibrium burnup. 
Note actual refueling upsurge to prepare for winter 1965— 
1966 peak-demand period. 


fuel were made in the following order: (1) fuel was 
moved only when the reactor needed reactivity, (2) 
channels with the highest burnup were refueled first, 
(3) problems unrelated directly to the fuel were con- 
sidered, and (4) flux profile was checked. The flux 
profile did not cause much difficulty because the con- 
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servativeness of the reactor design allowed sufficient 
maneuvering room. 

There are not yet quantitative data to explain the 
difference between the calculated and the actual 
schedules in Fig. 7; however, the difference is caused 
by (1) optimistic refueling rules in the calculated 
case; (2) imperfect refueling-maching operation [re- 
sulting in increased use of enriched booster (93% en- 
riched UA1), more shutdown fueling than desired, 
premature discharge of fuel, and some preferential 
fueling of the core center to maximize reactivity 
gain]; (3) the experimental program resulting in 
more discharge than demanded by the schedule; (4) 
mechanical difficulties requiring that specific chan- 
nels be left untouched; and (5) one accidental addi- 
tion of light water—approximately 1000 gal—that 
reduced moderator purity from 99.8 to 99.1% D,O. 
Continued refueling-machine development may elim- 
inate most of these difficulties. 


FOR THE FUTURE 


Experience with NPD indicates that, during the 
period before equilibrium, on-power refueling of a 
large power reactor must be scheduled with the aid 
of a sophisticated computer program that calculates 
the flux profile and the burnup in individual chan- 
nels. But after equilibrium is achieved, much of the 
computer work can be eliminated because the flux 
profile will be relatively static. Then refueling can 
be decided on the basis of (1) bundle burnups calcu- 
lated with static flux profiles or (2) decreased coolant 
temperature at a pressure-tube outlet as an indica- 
tion of burnup of fuel bundles in that tube. 


* * * 


This article has been prepared with the coopera- 
tion and assistance of Messrs. McConnell, Fanjoy, 
and Elston of the Hydro-Electric Power Commission 
of Ontario. 
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ENERGY CONVERSION 


Recent Developments in Liquid-Metal 
MHD Power Systems 


By Michael Petrick 


The growing interest in magnetohydrodynamic (MHD) 
generation of electric power stems from the potential 
15 to 25% increase in power-station overall effi- 
ciency possible with an MHD generator in a topping 
cycle on a steam plant.’ However, the modest tem- 
peratures achievable with economic central-station 
reactors likely to be available in the relatively near 
future make MHD generators of earlier interest for 
space applications. Recent developments in magneto- 
hydrodynamic technology, as described in the more 
than 200 papers presented in the third International 
Symposium on Magnetohydrodynamic Electrical Power 
Generation,’ focus attention on defining both the 
promise and the specific problems of the various 
types of MHD systems. This article summarizes re- 
cent significant advances in liquid-metal MHD cycles, 
which were discussed in detail for the first time at 
the third symposium. 

Future articles will similarly review other closed- 
cycle MHD systems and open-cycle MHD systems. 
Closed-cycle developments involved system experi- 
ments, nonequilibrium theory and measurements, con- 
formal mapping of the potential field and electrical 
leakage between electrodes. Discussions of open-cycle 
MHD concentrated on magnet and materials problems, 
striated flows, and combustion plasma physics. Gen- 
eral MHD problems, such as diagnostics, shock tubes, 
flow and field distributions, auxiliary ionization, and 
relaxation, were also discussed at the symposium. 


ROLE OF REACTORS 


The future of MHD power generation is very closely 
related to the direction and rate of development of 
nuclear power reactors. 

Closed-cycle MHD systems can become successful 
only if high-temperature reactors or equivalent heat 
sources can be developed. Reactors for closed-cycle 
gas plasma systems must operate at ~3000 °F, whereas 
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liquid-metal cycles require reactor temperatures 
1000 to 1500°F less. For either system the tempera- 
ture imposes specialized requirements on reactor 
design and development. 

On the other hand, and perhaps most important, the 
economic justification for MHD will be affected by 
the impact nuclear power is having on power costs. 
If nuclear fuel costs continue to decline as they are 
now doing, then the incentive for improving power- 
plant efficiency by means of MHD or other topping 
cycles will diminish. 


Status of Liquid-Metal MHD 


The feasibility of liquid-metal MHD conversion 
systems has been demonstrated. On the basis of the 
indicated performance of known fluid-flow devices 
(nozzles, separators, etc.), practical cycle efficiency 
of ~6% can be achieved. Such efficiency appears ade- 
quate for space power—and the general consensus 
is that liquid-metal MHD will probably be used ini- 
tially in various space applications. In fact, liquid- 
metal MHD power cycles could be available for use 
in space before turboelectric systems. 

There is considerable doubt, however, concerning 
the adaptability of liquid-metal MHD power cycles to 
commercial power systems. Minimum MHD topping- 
cycle efficiencies of 12 to 15% will be required to 
gain an economic incentive. It is not clear that such 
efficiencies can, in fact, be realized because of ex- 
tensive friction losses in various system compo- 
nents —especially in the 1100 to 1600°F range, which 
appears to be a reasonable extrapolation of today’s 
state of the art in reactor and materials technology. 

The most crucial area where additional work is 
necessary in liquid-metal MHD is in the development 
of efficient liquid accelerators for converting the 
thermal energy of the fluid emerging from the heat 
source to kinetic energy (or stagnation head) of the 








188 POWER REACTOR TECHNOLOGY 


fluid before it enters the MHD generator. The basic 
MHD cycles differ essentially in the method chosen 
for accomplishing this transformation. Generator de- 
velopment has proceeded more rapidly and is closer 
to resolution. There is little fundamentally new or 
different in liquid-metal MHD generators; the prob- 
lem is to tailor the generator to the cycle and achieve 
acceptable performance. As with the cycle, generator 
performance is also friction limited. A complete 
small-scale liquid-metal MHD power cycle could be 
operating at temperature within the next 5 years. 

Progress in liquid-metal MHD during the past 
2 years has been very substantial. The symposium 
called attention to the fact that liquid-metal MHD has 
emerged from its embryonic stage of development 
and is now a bona fide area of magnetohydrodynamics 
which is sparking considerable and increasing in- 
terest. 

The major areas of liquid-metal MHD which have 
received and are continuing to receive the most in- 
tensive study are: (1) generator performance, (2) 
methods for liquid acceleration, and (3) cycles. 


Generators 


Impressive progress has been made toward the de- 
velopment of liquid-metal MHD generators. Most of 
this has concerned either the induction generator 
operating with single-phase flow or the conduction 
generator operating with two-phase flow. 


A-C INDUCTION GENERATORS 

Progress has perhaps been greatest with induc- 
tion-generator development. An induction generator, 
such as is shown in Fig. 1, is a linear version of the 
conventional induction machine.* A traveling magnetic 
field is produced by polyphase windings, and, for the 
machine to run as a generator, the fluid velocity v 
must exceed the wave velocity v;, or the slip, defined 
as s = (v, —v)/vs, must be negative. Two phases are 
shown in Fig. 1, but a three-phase winding would 
probably be used, as in conventional machines. 


Sequence. Significant events in the development of 
induction generators occurred in the following se- 
quence. In the spring of 1964 the first induction gen- 
erator was tested at the Massachusetts Institute of 
Technology (MIT) by Jackson, Pierson, and Porter.® 
In August of 1965 a 350-cps 10-pole-pair flat-channel 
constant-fluid-velocity induction generator was tested 
by Prem of Atomics International on a NaK blowdown 
facility at Jet Propulsion Laboratory." Net power of 
1840 watts at 215 volts was generated with an electri- 
cal efficiency of 36%. Because operation was on a 
blowdown facility, the experiment ran for less than a 
minute. Evaluation of these tests shows that a sinus- 
oidal magnetic field and induced current and voltage 
can be produced. Extrapolation of these data to a 
large-scale generator led the author to conclude that 
efficiencies as high as 70% could be obtained. 
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The next induction generator was run in early 1966 
by Radebold and coworkers® at the Technical Uni- 
versity of Berlin. This generator produced a net 
power of 1 kw at 50 cps and 210 volts; the excitation 
power required was 9.5 kva. This generator, which 
ran for wl, hr under steady-state conditions, was de- 
signed primarily to produce basic data that could be 
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Fig. 1 Linear a-c induction MHD generator.° 


compared to theory, not to demonstrate high effi- 
ciency. 

Immediately preceding the symposium another sig- 
nificant development occurred which was not reported 
in a Symposium paper. The author had the opportunity 
of witnessing the startup of a generator from a ca- 
pacitor bank at the Technical University of Berlin. 
This successful experiment was also carried out by 
the AEG (Allgemeine Elektricitaits-Gesellschaft) group 
in collaboration with W. D. Jackson. 

Perhaps the most significant development in regard 
to induction generators was reported by Elliott and 
coworkers’—the development of a short (1 to 2) 
wavelength machine with compensation poles at the 
generator entrance and exit to reduce end losses. 
Generator development in the past progressed toward 
long-, many-wavelength machines to eliminate end 
effects. However, performance suffered because of 
extreme friction losses. The solution to the problem 
as proposed in Ref. 9 is to try to make an electrically 
short generator appear like a 1-wavelength section 
of an infinite generator. Analyses indicate that, with 
compensation poles, efficiencies up to 70% are at- 
tainable if the generator has a length/gap ratio of no 
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more than 70 (for low friction loss) and an electrical 
length of no more than 1 or 2 wavelengths (for low 
winding loss). A 1-wavelength generator has been 
built (Fig. 2) and is being tested at Jet Propulsion 
Laboratory.® 


Increasing Current Density. A way is being in- 
vestigated to overcome the potentially serious prob- 
lem of generator performance being limited by the 
current density in the generator windings.!® In MHD- 
generator configurations the windings system con- 
sists of a conventional multiphase system like that 
used for the conventional asynchronous machine. 
But, unlike these, the gap between iron cores is in- 
creased to a channel height generally 20 to 30 times 
greater. Therefore the reactance of the machine and 
the inductive component of the current become large, 
which in turn give the machine a small power factor. 
Compensation is made by an increase in the current 
density of the primary windings to about 20 to 30 
times that in the conventional asynchronous machine. 
This is only possible by use of a cryogenically cooled 
primary winding system. An MHD converter utilizing 
windings cryogenically cooled with liquid nitrogen and 
employing mercury as the working fluid is now being 
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run to study the predicted performance and the tech- 
nological problems involved. 


Materials Considerations. The generator channel 
must be able to withstand >1100 to 1600°F liquid 
metals flowing at several hundred feet per second 
with negligible corrosion or erosion." Refractory 
metals such as niobium and its alloys, which are 
good electrical and thermal conductors, are envi- 
sioned as candidates. The protection of the magnetic 
core and windings from excessive temperatures would 
require thermal insulation between the hot liquid 
metal and the field structure. Even so, to minimize 
insulation thickness and refrigeration power required 
for space applications, the magnetic core should be 
operated at the highest temperature consistent with 
high permeability and low loss. Thus the authors dis- 
cussed feasible operating temperatures for magnetic 
and excitation-winding materials, as well as cooling 
methods for the field structure. 

In considering design-engineering aspects, the au- 
thors extended MHD induction-generator theory to 
include both the conducting channel walls and the 
thermal insulation. Although the thermal insulation 
would have only a small effect on machine electrical 





Fig. 2 Experimental 1-wavelength generator is being tested at Jet Propulsion Laboratory. 
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performance, the channel wall thickness required (on 
the basis of pressure-vessel codes) to contain pres- 
sures of several hundred pounds per square inch 
would cause excessive electrical losses if conven- 
tional construction methods were employed. Thus 
novel construction schemes are suggested along with 
practical materials questions. (Some cycles do not 
involve generator-channel pressures requiring thick 
walls; for instance, the concept being pursued at 
Atomics International has 1- to 5-psi channel pres- 
sure.) 


D-C CONDUCTION GENERATORS 


The possibility of operating conduction generators 
on two-phase mixtures was discussed in Refs. 12 
and 13. Initial performance data on both a film-flow 
and a two-phase-flow conduction generator indicate 
that there are limited parameter areas of operation 
which are feasible.” As shown in Fig. 3, data ob- 
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Fig. 3 Efficiency decreases ina two-phase-flow genera- 
tor as mixture quality and void-volume fraction increase. 


Magnetic field was 4900 to 5500 gauss.4° 


tained from a two-phase (NaK—N,)-flow generator 
indicate that, as the mixture quality and void-volume 
fraction are increased, generator efficiency de- 
creases. The efficiency decrease is gradual until a 
mixture void fraction of ~0.85 is reached, then power 
and efficiency suddenly drop to virtually zero. The 
performance drop-off was attributed to rapidly de- 
creasing electrical conductivity of the two-phase 
mixture with increasing void fraction. 

The problem of defining the conductivity of two- 
phase mixtures remains unsettled, as evidenced by 
newly reported data for potassium mixtures at ele- 
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vated temperatures. Experiments carried out with 
potassium mixtures (vapor content varied from 2 to 
95%) at 1100°F indicated substantial electrical con- 
ductivity. The authors hypothesize that the emission 
of electrons from the surface of the liquid droplets — 
and not the influence of the liquid film —is the cause 
of the enhanced electrical conductivity. 

Preliminary film-flow-generator data were taken 
on a small-scale generator; the maximum measured 
power achieved was 237 watts. Performance was 
poor at very low liquid flow rates and magnetic fields 
but increased with increasing values of these param- 
eters. The limited data presented indicate that a 
film-flow generator apparently is feasible and stable 
in operation. However, many unanswered questions 
remain in regard to efficiency, range of stability, etc. 


REDUCING END LOSSES 


The end-loss problem has been analyzed for both 
a-c induction generators?,!4 and d-c conduction gen- 
erators.° End losses in a generator usually arise 
from the large gradients of magnetic field at the inlet 
and outlet of the generator channel. 

Methods were investigated for reducing the end 
losses in conduction generators, in both constant- and 
variable-area channels, by employing insulating vanes 
parallel to the flow.’® In an initial attempt to deter- 
mine optimum vane configuration, considering both 
electrical and viscous losses, a number of channel 
vane configurations were studied. The author was 
also able to extend the method of Sutton, Hurwitz, 
and Poritsky'® to channels of linearly variable area 
with the assumption of large Reynolds number and 
small magnetic Reynolds number. Under those con- 
ditions, and with the use of an “equivalent channel 
geometry,” a conformal-mapping technique permits 
the exact calculation of the net electrode current. 
For a 2:1 aspect-ratio channel with which no attempt 
has been made to use vanes to reduce end losses, 
maximum possible generator efficiency is 40.4%. If 
simple vanes are used, they must be very long, which 
means large viscous losses. Thus, with long vanes 
down the center line of the previous duct (both up- 
stream and downstream so that the aspect ratio be- 
comes 4:1), top efficiency rises to only 52,1%. Short 
but more complex vanes could yield 66% efficiency 
in the most promising configuration with a 2:1 as- 
pect ratio or 74% in a 4:1 channel. 

The scheme? proposed for resolving the a-c induc- 
tion generator end-effect problem by using compen- 
sating poles at the generator entrance and exit was 
discussed in the section on induction-generator de- 
velopments. 

In addition, previous analyses!’.18 of end losses in 
linear induction generators have been extended to in- 
clude circle diagrams depicting real and reactive 
currents.’ The two-dimensional induction MHD gen- 
erator of finite length is treated as a boundary -value 
problem by a two-sided Laplace transform. 
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MHD Cycles and Acceleration Methods 


Because of detailed analyses by numerous investi- 
gators, there is more general agreement on the pre- 
dicted efficiencies of the various proposed cycles 
than has ever existed in the past. All the papers on 
this topic dealt with three basic cycles, namely, the 
separator cycle, condensing injector cycle, and the 
two-phase generator cycle. The first two have been 
discussed and described in detail in the literature, '° 
whereas the latter was discussed in detail for the 
first time at this symposium.?).21 The two-phase gen- 
erator cycle was discussed initially in a preliminary 
manner at a previous symposium” and in Ref. 22. 


TWO-PHASE GENERATOR CYCLE 


This cycle can operate with either one- or two- 
component fluids (Fig. 4). In the one-component fluid 
cycle, the two-phase mixture generated in the heat 
source passes directly through the MHD generator to 
a condenser —diffuser and back to the heat source. In 
the two-component flow cycle, additional components 
needed are the mixer and separator. 

If the fluid is expanded through the MHD generator 
directly, the separate nozzle or stagnation head gen- 
erator can be eliminated. In a sense, the MHD gen- 
erator behaves as an expansion engine and produces 
electrical energy as the expanding conductive fluid 
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Fig. 4 Two-phase generator cycle variations can be either 
two-component system (above) or one-component system 
(below). 
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interacts with the magnetic field. The two-phase flow 
generator behaves in a manner analogous to that of 
the plasma flow generator. 

The efficiency of a two-phase generator cycle is 
potentially higher because power can be extracted 
under controlled velocity conditions to minimize 
friction—and the two-component mixture then can 
be separated at a lower energy and with lower losses. 
The potentially attractive performance of the cycle 
was discussed by Hammitt.”” He showed that higher 
cycle efficiencies are obtainable if a high gas-to- 
liquid volume flow ratio can be tolerated in the 
generator without undue loss in efficiency; under this 
condition, cycles that can use successfully a high 
volume flow ratio appear best. However, if the gen- 
erator efficiency decreases seriously with volume 
flow ratios, this conclusion is reversed. 

A novel two-phase generator cycle was also pro- 
posed by Bidard and Sterlini.”! The cycle features an 
MHD pump instead of a diffuser to return the fluid 
to the high-pressure end of the cycle. Another unique 
feature of this proposed cycle is that the entire ap- 
paratus would be enclosed within a magnet so that 
end effects in the pump and generator could be 
eliminated. 


CONDENSING-INJECTOR CYCLE 


Several new variations of the condensing-injector 
cycle were also proposed. Both Prem’ and Radebold® 
stress the possibility of minimizing the losses asso- 
ciated with mixing fluids of different velocities by the 
use of a multistage concept in which the injected- 
fluid velocity would be matched stepwise with the flow 
velocity in the nozzle. The theoretical thermodynamic 
efficiency of such a cycle was reported to be 60% of 
the Carnot efficiency,’ although no numerical analysis 
was presented to corroborate quantitatively the im- 
proved performance claimed. Present and planned 
studies will determine the conversion of a liquid into 
a foglike two-phase flow in the heat source, disper- 
sion of the liquid injected in the different steps of the 
transformation, and condensation of the fog flow. 

Because of the potential for considerably higher 
internal efficiency, the two-component cycle was 
perhaps the most important variation of the condens- 
ing injector analyzed.*3.24 With lithium as the sub- 
cooled injection stream and cesium or potassium as 
the vaporized stream, the losses due to the accelera- 
tion of the injection stream would be reduced because 
of the high heat-sink capacity (and thus lower injec- 
tion flow rate) of lithium. For a cycle with lithium — 
cesium as the two-component working fluid, vapor — 
liquid internal efficiency may approach 80% and cycle 
thermal efficiency could reach 30% with a 1200°K 
maximum cycle temperature. There is an optimum 
vapor quality that gives maximum injector and cycle 
performance; lower mixing-chamber pressure re- 
duces injector efficiency but raises overall effi- 
ciency; increased inlet liquid velocity diminishes 
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performance. Because the results presented in both 
studies are based solely on thermodynamic consider- 
ations and no losses are accounted for, the results 
represent the maximum thermodynamic performance 
that can be achieved. Although the results are com- 
parable to those presented previously for other fluids 
and cycles,?5—27 these studies underscore the possible 
gains that might be obtained through the use of two 
fluids in the cycle. 

Additional experimental data”® on the performance 
of the condensing injector for space power show exit 
pressures greater than both liquid and vapor inlet 
pressures. The data, which are encouraging in terms 
of predictability and efficiency performance, also 
indicate that detailed knowledge of the mixing pro- 
cess is needed to design high-performance condens- 
ing injectors. The highest efficiencies occur when 
the velocity ratio between vapor and liquid is maxi- 
mum and the flow ratio of liquid to vapor is mini- 
mized, 


Progress Expected in Next 2 Years 


The greatest effort in the next 2 years will probably 
be directed toward building and testing an efficient 
induction generator. As rapporteur, W. D. Jackson 
stated in his analysis of the state of the art of gen- 
erator development, “A quick calculation of the rele- 
vant magnetic Reynolds number and a comparison of 
it with what rotating machine designers call the 
reactance to the distance ratio allowable in solid 
rotor induction machines shows that the results pre- 
sented for the operation of these generators, their 
self-excitation, their steady running, and so on, was 
never really seriously in doubt. The basic issue re- 
mains one of realizing an induction generator of 
acceptable real efficiency within the limitations im- 
posed by systems and materials considerations.” 

Another major question that is expected to be 
resolved in this time period will be the feasibility 
and performance of the condensing-injector cycle. 
The geometries, fluids, and parameters of interest to 
liquid-metal power cycles have not been studied ex- 
perimentally thus far. However, studies under way at 
MIT and Argonne National Laboratory and elsewhere 
will remove the uncertainty concerning injector per- 
formance. 
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Component Development: 
Key to Rankine-Cycle Space Applications 


By Edward J. Croke 


The development of reactor and isotope heat sources, 
combined with projected requirements for compact, 
lightweight electric-power generators for space and 
undersea application, has given rise to extensive 
turboelectric-power development programs intended 
to exploit the potentially high efficiency of the Ran- 
kine cycle. Experience in these system development 
programs has spotlighted the importance of im- 
proving the technology of high-performance com- 
ponents capable of long-time unattended operation. 

In terms of practical engineering, a great deal of 
Rankine-cycle component development has taken place 
in recent years. Much of it has been based on very 
pragmatic developmental methods. At present, bear- 
ing wear, boiler conditioning, and turbine erosion 
are the most significant problems in achieving reli- 
able Rankine cycles. Fabrication techniques for 
longer life bearings and higher temperature turbines, 
more efficient pumps, and less vulnerable radiators 
are being developed. Heat-transfer properties of 
useful liquid metals are being investigated in sev- 
eral research facilities. Nevertheless, for the long 
term, the work done so far indicates that additional 
basic research is needed in the areas of heat trans- 
fer, materials, and structural mechanics. The situa- 
tion is illustrated by consideration of the status of 
the four major Rankine-cycle programs and the com- 
ponent experience discussed in this article. 


Rankine-System Development Programs 


The four major Rankine-system development pro- 
grams include the mercury Rankine program spon- 
sored by the U. S. Atomic Energy Commission 
(AEC),! the SNAP-8 Program of the AEC and the 
National Aeronautics and Space Administration 
(NASA),? the SNAP-50/SPUR Program of the AEC 
and the U. S. Air Force,’ and the Organic Working 
Fluid Rankine Cycle Programs of the U. S. Navy, 
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Air Force, and AEC.‘ There also have been asso- 
ciated component studies and development efforts at 
various government and industrial facilities. 

The early history of Rankine-system development 
parallels the growth of the national space program. 
Initially, mission requirements were not well defined, 
and, with a premature sense of optimism, full-scale 
systems were projected before the problems of 
building a component technology were fully evident. 
When experience spotlighted the component problems, 
and aS mission requirements were seen to change 
rapidly, most of the systems programs evolved into 
component and technology programs that will permit 
development of a range of Rankine systems that can 
supply power for a variety of anticipated space mis- 
sions. 


MERCURY RANKINE PROGRAM 


Development of the AEC’s SNAP-2 system began 
in 1957 with the objective of producing a flyable 
3-kw(e) reactor-powered mercury Rankine system 
(Fig. 1) capable of operating for 1 year (Ref. 5). 
With changes in mission requirements, the SNAP-2 
program evolved into the mercury Rankine pro- 
gram-—whose initial objectives were to develop 
system technology and qualify components (exclusive 
of the reactor or isotope heat source) for a space 
power plant with 3- to 5-kw(e) power-conversion 
modules so that the plant could be adapted to supply 
3 to 25 kw(e) for 2500 hr of operation in space.! The 
ultimate goal was to achieve reliable performance 
for 10,000 hr. The mercury Rankine program is now 
in the process of being phased out, with only two 
combined-rotating-unit turbomachinery tests to be 
completed. 


Component Tests. Development tests of the com- 
ponents at Atomics International, the prime con- 
tractor, indicated that the required operating life 
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could be achieved. System performance, orbital start- 
up capability, and structural integrity have been 
verified in system-test facilities.' Although one com- 
bined rotating unit was tested for 7000 hr without 
major damage, the major developmental problem 
has been journal-bearing erosion, which could limit 
the life of the mercury rotating machinery. 

The previously experienced bearing erosion— 
cavitation has been decreased by almost an order of 
magnitude by reducing loads and increasing bearing 
pressure. Life capabilities well beyond 10,000 hr 
now seem feasible. 

Other components demonstrated satisfactory per- 
formance. The model-5C boiler, which was tested 
for 38,000 hr, has proved to be structurally and 
thermally acceptable for long-life operation if pro- 
tected from exposure to air and contaminants such 
as trace quantities of oil. The adequacy of the 
thermoelectric pump was demonstrated by many 
hours of testing. Several pumps have operated for 
>20,000 hr. Both the honeycomb and aluminum-fin 
radiator —condensers have demonstrated good struc- 
tural and operational characteristics, although the 
process for bonding aluminum radiator fins to the 
Haynes 25 condenser tubes caused some develop- 
mental problems. The startup components (pressure 
regulator and injection tanks) posed no problems 
more serious than the fabrication of the metal start- 
ing bellows. 


System Startup. Orbital (zero-gravity) transient 
characteristics of the integrated Rankine system 


1240°F 265 psia 
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have been simulated successfully, and automatic 
startup has been accomplished with a nonnuclear 
installation. Although the simulated system achieved 
stable operation, long operating lifetimes remain to 
be demonstrated —and performance degradation, con- 
trol interactions, mercury corrosion, metallurgical- 
processing effectiveness, and long-term exposure to 
vacuum must be evaluated to qualify the system for 
flight operation. Multiple-restart capability has also 
been demonstrated. 


SNAP-8 SYSTEM 


The NASA Lewis Research Center has project- 
management responsibility for the development of 
SNAP-8 power-conversion equipment, and the AEC 
Space Nuclear Systems Division is charged with 
nuclear-systems development. 

The mercury Rankine turboelectric-power system 
(Fig. 2) for the SNAP-8 Program, which was based 
considerably on technology developed under the 
SNAP-2 program, has had a history somewhat simi- 
lar to that of the SNAP-2 Program and the mercury 
Rankine program. The initial objective of developing 
a 35-kw(e) reactor-powered system for 10,000 hr of 
operation in space after a remote-controlled startup 
has remained more or less intact from program in- 
ception in 1960 to the present. However, limitations 
of the state of the art of liquid-metal-systems tech- 
nology in 1962 resulted in a major redesign of the 
prototype system. “With system-level development 
curtailment in 1965, actual development was limited 
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Fig. 2 The SNAP-8 mercury Rankine-cycle system, developed for NASA Lewis Research Center by Aero- 
jet-General Corporation, consists of a compact uranium—zirconium fueled-and-moderated reactor (devel- 
oped by Atomics International), a boiler, condenser, turboalternator, radiator—condenser, mercury and 
Nak pumps, and a control system. These are arranged in three liquid-metal circulating loops: A NaK pri- 
mary loop cools the reactor and powers the mercury Rankine-cycle secondary loop; the tertiary Nak loop 
transfers heat from the mercury condenser to the radiator—condenser, where the heat is rejected to space. 
In a fourth loop an organic fluid circulates to lubricate bearings and to cool the alternator, electric motors, 
and controls.” 
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to an attempt to qualify the first power-conversion 
test system for 10,000 hr of operation. 


Reactor Tests. The SNAP-8 test reactor went crit- 
ical in June 1963 and later demonstrated sustained 
power operation and static and dynamic stability. 
The longest continuous test run was 5000 hr; the 
reactor and test facility were shut down as planned 
in April 1965 after accumulating a total of 12,000 
test hours that included 100 days of operation at 
600 kw(t) with a NaK exit temperature of 1300°F. 
Design has been completed, and fabrication is under 
way for the second SNAP-8DR reactor test. 


Component Tests. Although the boiler designed for 
SNAP-8 proved generally satisfactory, it was plagued 
with inconsistent boiler conditioning (scouring of 
scale oxides and contaminants from the heat-transfer 
surface). The mercury condenser performed gen- 
erally well in preliminary tests. The 51.6-kw(e) 
turboalternator performed satisfactorily in more than 
830 hr of tests at approximately design conditions. 
The NaK pump-motor assemblies accumulated more 
than 1000 hr of operation in > 1100°F Nak; the mer- 
cury pump-motor assemblies completed > 120 hr of 
operation at an inlet temperature of 500°F; a Thomp- 
son Ramo Wooldridge Inc. lubricant—coolant pump- 
motor assembly achieved >1000 hr of operation at 
rated conditions without performance degradation 
(seven of these units accumulated a total of more 
than 4300 hr of operation). The start system, voltage 
regulator, speed controller, and related electrical 
subsystems demonstrated satisfactory performance. 

Although the mercury boiler, turbine, condenser, 
and pump demonstrated satisfactory combined per- 
formance for 837 hr in a multiple-component test 
loop, the entire nonnuclear four-loop system has not 
yet been operated as an integral unit. Such a test of 
the first power-conversion system is scheduled for 
early 1967 as part of the component qualification 
program. 


SNAP-50/SPUR SYSTEM 


Little information is available on this program, 
but AiResearch Mfg. Co. (a division of the Garrett 
Corp.) indicates that turbine-rotor bearing materi- 
als, lubrication, and stability may be among the 
foremost problems. °® Blade erosion and associated 
turbine-life limitations are no longer expected to 
loom large among development difficulties. Experi- 
ence with these components will be considered in 
greater detail in the Rankine Component Technology 
subsection of this article. 


ORGANIC RANKINE-CYCLE SYSTEM 

This 1.5-kw(e) Rankine system, whose combined 
rotating unit was being developed by Sundstrand 
Corp. for the U. S. Air Force, would have been 
heated either by a *!°Po source or by solar energy 
concentrated by a mirror. The Dowtherm A working 


fluid reaches 4 maximum temperature of ~700°F, 
and so the high-temperature problems of the liquid- 
metal Rankine systems are avoided. The organic 
Rankine system (Fig. 3) consists of a boiler, a com- 
bined rotating unit (including turboalternator and 
Pitot pump), a radiator—condenser, a regenerator, 
and the control system. The thermodynamic cycle 
is shown in Fig. 4. 

The combined rotating unit accumulated more than 
2060 hr of successful tests; a single set of bearings 
has served for more than 600 hr of cumulative run- 
ning time with more than 100 restarts. Although these 
tests have shown good performance, long-endurance 
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Fig. 3 Organic Rankine-cycle turboelectric system for 
Apollo type mission would have a 17.5% overall cycle ef- 
ficiency (including losses and control power) with 750°F 
157-psia Dowtherm A working fluid and *!°Po as fuel. The 
system would provide 1.5 kw(e) for a 14-day mission plus 
a 30-day hold. Fuel block and shielding would weigh 250 lb, 
and redundant conversion equipment would weigh 280 lb 
(Ref. 4). 
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Fig. 4 Temperalure—entropy diagram for Dowtherm A 
Rankine cycle. Dowtherm A has a composition of 26.5% by 
weight (CgH5). and 73.5% (CgHs)0, molecular weight of 
165.635, freezing point of 33.6°F, and vapor pressure of 
0.03 psia al 300°F and 106.81 at 700°F, Corresponding 
characteristics of biphenyl, another organic fluid, are 
(CsH5)2 composition, 154.20 molecular weight, 198°F freez- 
ing point, and vapor pressure of 0.68 psia at 300°F and 
109.0 psia at 700°F (Ref. 4). 
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system capability equivalent to 10,000 hr of unat- 
tended operation in space remains to be demon- 
strated. Because the thermal stability of the working 
fluid may limit system life, AEC programs are 
evaluating this uncertainty. 


Rankine Component Technology 


Each major component of a high-temperature Ran- 
kine system has its own developmental and opera- 
tional problems. Thus it is worthwhile to consider 
and compare components developed in the various 
Rankine programs. In general, parallel and essen- 
tially independent solutions seem characteristic of 
Rankine technology. 


PUMPS 


Notable among the pumps under development are 
the mercury and NaK pumps for SNAP-8, the NaK 
pumps for the mercury Rankine program, the potas- 
sium pumps operated at Oak Ridge National Labo- 
ratory (ORNL), and those being developed by General 
Electric Company for NASA. 


Mercury Pumps for SNAP-8. The duty cycle of the 
mercury pump in SNAP-8 includes 10,000 hr of un- 
attended operation in a radiation environment of 10'° 
neutrons/cm? at a speed of 8000 rpm and with an 
inlet suction temperature of 405°F. Pump flow 
rate is 1.775 gal/min at a nominal head of 10.5 psi 
(Ref. 2). 

A jet pump at the inlet of the main pump housing 
overcomes potential cavitation difficulties at the 
rather low specified pump suction head. The mercury 
pump is a single-shaft overhung system with an im- 
peller at one extremity (upstream of the jet pump) 
and a drive motor at the other. The semiopen-vane 
centrifugal impeller is ball-bearing mounted in a 
housing that has been designed to minimize bearing 
loads. The bearing lubricant, organic 4P3E (poly- 
phenyl ether), is separated from the mercury flow 
by a dynamic seal system consisting of a viscopump 
that maintains a stable liquid—vapor interface in the 
fluid. Leakage of the fluid to space is minimized by 
a helical thread molecular pump that has a 3-mil 
rotating clearance. 

On the bearing side of the rotating system, a 
Similar dynamic seal is used; however, additional 
low-speed (6000 rpm) and stationary sealing is pro- 
vided by stationary carbon face seals; a bellows 
arrangement lifts the carbon face seals off the bear- 
ing surface at rotational speeds > 6000 rpm. The 
central housing of the pump supports two angular 
contact end bearings of M-50 tool steel. Between 
these end bearings is a rotating “slinger” seal and a 
molecular pump. The pump motor has silicon-steel 
laminae and is a 400-cycle induction unit cooled with 
4P3E. At the ends of the motor, the terminals are 
hermetically sealed. Between the motor and one end 
bearing is a scavenging seal; at the other end of the 
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pump, there is a face seal. The prospect of running 
the pump flooded was considered, but the system 
temperature rise caused by viscous effects proved 
to be excessive, and hence the end scavenger was 
provided. 

The test pump weighs ~186 lb in the fully instru- 
mented state. Much of this weight is in the 9% chro- 
nium—1% molybdenum steel housings hogged from 
solid forgings. A realistic flight weight for the sys- 
tem would be ~110 lb. 

Tests have been conducted on both the subas- 
sembly and system levels. The cavitation perfor- 
mance of the jet pump has been tested with mercury, 
and the sealing system has been tested for 1000 hr 
in vacuum. Extrapolation of the results of these 
tests shows that the mercury loss to space will be 
~1.75 lb per 10,000 hr and that 0.1 lb of 4P3E will 
be lost in this same time. The motor has been tested 
for over 20,000 hr immersed in 4P3E, and, although 
the coolant temperature reached 350°F during the 
tests, no compatibility problems were evident. 

Three pump assemblies were built. Two were 
retained at Aerojet-General Nucleonics and one was 
sent to NASA Lewis Research Center for testing. 
More than 140 starts were accumulated on these 
units, of which 72 starts were on one pump. Except 
for two failures of the space-seal lift-off bellows, 
all tests were successful. However, some degree of 
cavitation was detected near the end of the tests. 

The cavitation limit of the assembly was imposed 
by the centrifugal impeller, which is provided with 
a 2-ft head by the jet pump. The minimum net posi- 
tive suction head ranged from 0.1 ft at 0.5 gal/min 
to 0.75 ft at 2.25 gal/min. Overall efficiency of the 
pump was found to be ~13%. 


Nak Pumps for SNAP-8. In the primary and ter- 
tiary NaK loops of the SNAP-8 system, each pump 
is a self-contained hermetically sealed unit with a 
single-shaft connecting motor (208-volt induction, 
rated at 1530 rpm) and an impeller. The type 410 
stainless-steel impeller, a semiopen-vane unit with 
90° backvanes, has a rated specific speed of 
1600 rpm and is enclosed in a double volute housing 
of type 304 stainless steel. It is designed for use 
with 1170°F fluid.’ 

The pump has been tested in a loop containing a 
heat-exchanger cold trap, an economizer, and a 5-u 
NaK filter. The fluid discharged at 350°F into the 
economizer, where it was cooled to 320°F and sent 
to the heat-exchanger cold trap and cooled to 290°F. 
It then entered a canister containing stainless-steel 
wire mesh (to trap oxides) and returned to the econ- 
omizer, where it was heated to 300°F and filtered. 
(Reheating prevents the filter from becoming an oxide 
trap since it is not permitted to become the coldest 
part of the system.) 

Considerable testing of this design also has been 
done on the subassembly level. The motor has been 
dynamometer checked in air, and the bearings have 
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been run in both oil and NaK to determine the mini- 
mum allowable bearing thickness as a function of 
bearing load. The capability of the cold trap to keep 
the oxide level below 20 ppm was tested and verified, 
and water tests have determined the heat-exchanger 
heat-transfer characteristics. 

A major development problem was NaK leakage 
(at temperatures above 500°F) through inclusions in 
the original type 347 stainless-steel pump housing. 
This design was successfully fabricated in type 304 
stainless steel. Another problem was rotor sticking, 
caused by flat bearing pads; this situation was cor- 
rected by grinding a radial taper on the bearings to 
obtain a line contact. 

Eight pumps were built and tested at fluid tem- 
peratures up to 1200°F. One unit was run for 2040 hr. 
The design is reported to have met or exceeded all 
requirements with a demonstrated efficiency of 35 
to 40%. (The significant NaK pump losses are fric- 
tional.) 


Thermoelectric-Powered Nak Pump. The Atomics 
International electromagnetic pump (Fig. 5) is an 
unusual device in that the 3.5 kw(e) required for its 
operation is generated thermoelectrically from the 
waste heat of the primary loop.® Rated for a flow of 
12 gal/min at 1200°F, the pump has a pressure head 
of 1 psi. The design has been optimized in terms of 
weight vs. thermoelectric thickness, weight vs. throat 
height, and weight vs. magnetic-field strength. 


Pumped fluid Magnetic 
(hot body) field 
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Fig. 5 Schematic of SNAP-2 thermoelectric-powered elec- 
tromagnetic NaK pump.’ 


An electrically heated test loop was constructed 
for this direct-current induction pump. Four test 
series were devised: (1) flight-operation duty cycle, 
(2) startup transient tests, (3) steady-state perfor- 
mance evaluation, and (4) endurance testing. Of the 
12 pumps built, 10 have been tested at rated condi- 
tions. About 20,000 hr have been accumulated on each 
of three units; two tests were terminated after 
5000 hr of operation each. 

Tr: maximum temperature transient encountered 
was 363°F/sec—well above the specification limit 
of 50°F/sec —but no adverse results were reported. 
Maximum observed degradation in performance was 
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associated with a 16% flow decrease after 5000 hr of 
operation on one unit and a 10% flow decrease after 
5000 hr on another. Two units, tested in a vacuum 
environment for 800 to 1000 hr, showed no signifi- 
cant performance degradation; test units subjected 
to thermal cycling evidenced leakage only after 3000 
to 5000 hr. Magnet degradation caused by nuclear 
radiation at operating levels appeared negligible 
over the time considered. 


Potassium Pumps with Turbine Drives. A program 
was established at ORNL to accumulate experience 
with liquid-metal turbopumps by fabricating and 
testing several machines of typical performance and 
behavior.? High reliability was found to depend on 
bearing design and materials, minimization of tur- 
bine-bucket damage, and minimization of pump cav- 
itation. 

The turbopump developed at ORNL (Fig. 6) used 
bearings of tungsten carbide with 12% cobalt binder. 
Because the coefficient of thermal expansion of this 
material differs significantly from that of the aus- 
tenitic steel casing, TZM bearing housings were 
chosen. Radial thermal expansion was accommodated 
by cantilever-mounting the bearings to obtain radi- 
ally soft but axially stiff behavior. Journal bearings 
were lubricated with potassium. Axial loads were 
carried by two flat-plate thrust bearings. 

The turbine shown in Fig. 6 is 2.75 in. in diameter, 
and the impeller is 1.0 in. in diameter. Turbine 
speed is 24,000 rpm, and the design inlet tempera- 
ture is 1540°F. Eight of these units were fabricated, 
and seven were successfully tested to failure for 
periods ranging from 0.5 to 327 hr. One continuous 
run was 140 hr, and total accumulated operating 
time was >540 hr. More than 125 starts were re- 
corded. In these tests, termination was caused by 
such failures as bearing-sleeve fractures, seizing 
of the journal bearings of the turbine or impeller, 
impeller rubbing against the pump housing, or filter 
clogging. Figure 7 shows pump performance. 

Fabrication of two larger pumps was contracted to 
Aeronutronic/Philco. The first of these, Mark I 
(Fig. 8), had a 6-in.-diameter turbine wheel and 
a 2.75-in.-diameter impeller; the second model, 
Mark II, was similar but was rated for a higher flow 
and had a more massive casing. Other improve- 
ments in the Mark II included (1) carburizing all 
TZM surfaces to eliminate galling, (2) nitriding 
all stainless-steel surfaces to reduce corrosion, 
(3) mounting the WC bearing housings on conical 
surfaces to compensate partially for differential 
thermal expansion, and (4) lengthening certain bear- 
ings. The Mark I turbine used convergent —divergent 
nozzles, but the Mark II turbine used only con- 
vergent nozzles. 

Characteristics and performance of three turbo- 
pumps are shown in Table 1. Although one of these 
pumps was run to 75% of theoretical cavitation con- 
ditions, cavitation damage was not found. The low 
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Fig. 6 Turbine-driven potassium 
boiler feed pump.® 
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Fig. 7 Turbine-driven-pump performance characteristics 
on steam and water. Pr is turbine inlet pressure, W is 


turbine vapor flow; turbine outlet pressure is ~2 psia for 


all runs.9 
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Pump discharge 


velocities and short runs may have prevented dam- 
age, but there is reason to believe that extreme 
cavitation need be no more damaging to a machine 
than mild cavitation—depending upon the properties 
of the fluid and the degree of subcooling involved. 
The quality of the turbine exhaust was ~ 95% for these 
units. Little erosion was detected under these con- 
ditions, in which the condensated consisted of very 
finely divided particles. 


BOILERS 


A principal factor limiting Rankine-system boiler 
design has been a lack of precise liquid-metal heat- 
transfer data. A better design basis has been es- 
tablished by studies with mercury, potassium, and 
sodium. 


MRP Mercury Boiler. Considerable test experi- 
ence has been achieved with the mercury Rankine 
program compact boiler; three design models have 
accumulated a total of 38,000 hr. Major problems 
have been with the mercury corrosion and boiling 
heat-transfer conditioning of the Haynes 25 tubing. 
The mercury corrosion rate is diffusion controlled, 
decaying with time such that the generated “crud” 
products can be adequately trapped or filtered to as- 
sure at least 1 year of operation. The boiler decon- 
ditioning appears to be controllable if precautions 
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are taken with respect to cleanliness, air inclusion, 
and surface treatment. Even with some degree of 
boiler deconditioning, the power from the model-V 
combined rotating unit was generally adequate. The 
endurance experience of various boiler tests indicates 
generally that conditioning characteristics improve 
with time with a leaktight system. 


Mercury Boiler for SNAP-8. Considerable boiler 
test experience was accumulated during the SNAP-8 
development effort.'° A 480-kw(t) forced-convection 


Table 1 TURBINE-DRIVEN POTASSIUM BOILER FEED PUMPS? 





Pump 


ORNL Mark I Mark II 





Speed, rpm 25,000 8,000 8,400 
Inlet temp., °F 1,540 1,540 1,540 
Turbine inlet pressure, psi 30 30 30 
Turbine exhaust temp., °F 1,050 1,050 1,140 
Turbine exhaust pressure, psi 1.5 1.5 3.0 
Flow, gal/min 1.5 1.5 3.0 
Suction temp., °F 1,050 1,050 1,050 
Suction pressure, psi 1.5 1-4 1-4 
Pump discharge pressure, psi 70 50 60 
Accumulated run time, hr 327 1,500 








Pump discharge 


Mark I turbine pump had convergent—divergent nozzles.9 


one-pass mercury boiler heated by NaK (with a gas- 
fired heater) was extensively tested. The NaK enters 
this tube-in-shell boiler at 1280 to 1330°F. Mercury 
enters (0.3 ft/sec) the double spiral helical coils 
surrounding the NaK shell at 320 psia and 517°F, 
and, at rated operation, exits at 1252°F. The helical 
tubes are 9% chromium—1% molybdenum steel, chosen 
for its corrosion resistance, and the NaK shell is of 
type 315 stainless steel. Over 1600 test hours were 
accumulated on this configuration —1400 hr were on 
a single unit. 

These boilers are rated for 11,000 lb/hr flow and 
have mass velocities of 200 lb/(sec)(sq ft). Insta- 
bility was 4 to 6 cycles/min or ~ 200 psi. Predictions 
of pressure drops and heat transfer were made with 
the Martinelli correlation. 

Boiler conditioning proved to have a major in- 
fluence on performance. The first unit tested re- 
quired ~ 100 hr of operation to become wetted and to 
achieve rated performance. When conditioning was 
complete, the pressure drop and the exit tempera- 
ture stabilized near the predicted values for a wetted 
system. A second boiler operated for ~20 hr in a 
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derated condition with an exit temperature of ~ 900°F 
rather than the rated 1252°F. Exit quality was only 
~70%, and the pressure drop was correspondingly 
low. After 20 hr the exit temperature stabilized at 
1200°F, and near rated performance was achieved. 


Later tests revealed that the addition of ~20 ppm 
liquid rubidium (to serve as a wetting agent) rapidly 
raised the exit temperature of an unconditioned 
boiler from 800 to 1200°F with a corresponding in- 
crease of pressure drop. After boiler shutdown and 
restart, further rubidium injection was required — 
probably because of consumption as oxides during 
the shutdown period. Because rubidium does form 
oxides that cause unpredictable alteration of the 
heat-transfer properties of the mercury, this con- 
ditioning technique was not favored for the SNAP-8 
system. After 200 hr of operation with rubidium, the 
boiler was shut down, cleaned out, and restarted; 
rated performance was attained immediately. Ru- 
bidium appears to have promoted the cleaning of 
scale from the boiler-tube surfaces. 


Heal Transfer in Mercury Boilers. Four basic 
modes of heat transfer in mercury flow have been 
defined: (1) contact boiling in the liquid regime, (2) 
slug flow, (3) intermediate contact (semiboiling), 
and (4) film boiling. The first, third, and fourth of 
these flow regimes are fairly well understood, but 
the second, slug flow, is neither well known nor easy 
to investigate. Liquid slugs form at low flow veloci- 
ties (~0.5 ft/sec) and tend to use up considerable 
boiling length without significant heat transfer. High 
flow velocities are required to avoid slug flow ina 
heat exchanger; however, the associated pressure 
drops are frequently prohibitive. To achieve a high 
pressure drop and flow velocity in a test boiler at 
Aerojet-General Nucleonics, workers fabricated sev- 
eral inlet plugs that would provide high flow veloci- 
ties.'' These were constructed in two sections, a 
high-velocity section that gave 0.9 to 6 ft/sec with 
some preheating of the fluid and a two-phase section 
with a swirl inducer. Each plug consisted of a cylin- 
drical core spirally wrapped with wire. The amount 
of flow constriction was such as to initiate boiling 
at the exit of the first section of the plug since, if 
boiling were not in progress at the plug exit, the 
high velocities achieved would not be maintained 
downstream of the plug, and slug flow might ensue. 


A plug designed for an exit velocity of 0.9 ft/sec 
was installed in a helically wound NaK—mercury 
boiler fabricated from type 315 stainless-steel tubing. 
Tests with this plug indicated that the boiler re- 
mzined unconditioned over most of its length even 
after extended operation. With a 6 ft/sec plug, sig- 
nificant conditioning was observed after 4 hr and 
nearly complete conditioning was observed after 
20 hr of operation. Subsequent tests with a 4.5 ft/sec 
plug and with a 0.9 ft/sec plug indicated full condi- 
tioning. From these tests was drawn the rather ob- 
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vious conclusion that a fully conditioned boiler is 
free of slug flow and needs no inlet plug. 

Contaminated fluid tests run with the same set of 
plugs, but with the addition of oil (contaminant) to 
the fluid, proved more interesting. When the condi- 
tioned boiler with the 0.9 ft/sec plug was contami- 
nated with oil, no deconditioning was observed until 
the system was shut down. When restart was at- 
tempted, the boiler was found to be deconditioned; 
no amount of running would recondition it. 

A conditioned boiler was then fitted with the 
4.5 ft/sec plug. The addition of oil caused decondi- 
tioning, but reconditioning began after continued 
operation. Regression was observed after shutdown, 
but full reconditioning was achieved after restart 
and 75 hr of operation. In the next test the boiler 
was fitted with the 6.5 ft/sec plug, and oil was added. 
Deconditioning occurred, but reconditioning was com- 
plete in only 25 hr. Higher plug velocities promoted 
corrosion and erosion in the systems. (Each of the 
test plugs was designed for 30% exit quality.) 

Deconditioning inevitably occurred after boiler 
shutdown in these tests. Thompson Ramo Wooldridge 
reported that hydrogen bright annealing of boilers 
tends to minimize or to eliminate this effect. 


Potassium Boiler for SNAP-50/SPUR. A multitube 
potassium boiler, designed to test overall heat- 
transfer performance against theoretical predictions 
and to verify the results of Freon tests, was fabri- 
cated at AiResearch Mfg. Co.* The unit, a 19-tube 
Haynes 25 tube-in-shell boiler, is counterflow heated 
by sodium in the primary loop. The outer shell is 
semicircular; the inner 0.25-in.-OD tubes contain 
twisted tape inserts diffusion bonded to the tube walls. 
Potassium is circulated by an induction pump; an 
auxiliary heater at the boiler inlet assures the de- 
sired fluid inlet temperature. 

The design flow was 5 to 15 lb/min at a mass 
velocity of 10 to 30 lb/(sec)(sq ft). Exit conditions 
of up to 150°F superheat were anticipated. Tests 
indicated that system performance was correlated 
well by means of the Martinelli equation incre- 
mented by 30% to correct for the twisted tape in- 
serts. For fluid exit qualities as low as 40%, the 
overall measured heat-transfer coefficient was lower 
than the predicted value, as for low flow rates 
(~5 lb/sec), where data scatter was considerable. 

Operation of the boiler was stable at all times, 


‘and boiling was easily initiated by providing sodium 


at only ~10°F above the potassium saturation tem- 
perature. Superheating up to 150°F was readily at- 
tained. The achievement of a pressure drop sufficient 
to ensure two-phase flow was achieved with inlet 
orificing. 


High-Performance Potassium Boiling. A 300-kw 
Haynes 25 potassium boiler was fabricated by Gen- 
eral Electric Company in Evendale, Ohio. This unit 
employs a gas-fired sodium loop as the heat source 
and features an air-cooled condenser. The test sec- 
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tion is illustrated in Fig. 9. A series of tests mea- 
sured local heat-transfer coefficients over the entire 
series of flow regimes from the subcooled liquid to 
the superheated vapor states. Local heat-transfer 
data were evaluated in the nucleate-boiling, critical- 
heat-flux, transition-boiling, and superheated vapor 
phases. The tests also produced average heat-trans- 
fer data and pressure-drop measurements.'” 

Data runs included tests in which subcooled heat- 
ing, nucleate boiling, transition-zone boiling, and 
superheating occurred simultaneously in the flow 
channel. Local heat fluxes and heat-transfer coef- 
ficients were measured in each zone. Figures 10 to 
12 are typical of the results obtained. This infor- 
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Fig. 9 Instrumentation of 300-kw test facility that includes 
Haynes 25 potassium boiler.” 
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Fig. 10 Results of a typical once-through boiling run with 
P/D = 6 insert in the 300-kw facility shown in Fig. 9. Data 
points show calculated or known values; solid lines shou 
average values; dashed lines show estimated values.” 
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mation should be useful in the design of potassium 
boilers over a wide range of operating conditions. 

Additional nucleate-boiling data were taken at a 
power level of 100 kw. The effect of vapor quality on 
heat transfer was not observed to be significant, but 
the magnitude of the total heat flux did influence the 
heat-transfer coefficient. 
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Fig. 11 Critical-heat-flux results for a 0.67-in.-ID boiler 
tube, showing (top) the effect of mass velocity and (bottom) 
the effect of temperature.” 
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Pressure-drop data for boiling potassium were 
taken at 300 kw (Fig. 13); fairly good agreement was 
established with the Martinelli correlation. It ap- 
peared that the variation of the overall heat-transfer 
coefficient could be predicted with reasonable ac- 
curacy for flow conditions which were not very dif- 
ferent from those set up for the tests, but, for flow 
conditions which differ significantly from those tested, 
local heat-transfer coefficients must be employed. 
A homogeneous flow model usually predicts a lower 
heat-transfer coefficient than is actually the case. 


Forced-F low Boiling of Sodium. A high-tempera- 
ture test loop was assembled at NASA Lewis Re- 
search Center to study the heat-transfer character- 
istics of boiling sodium.'* All components, including 
the 4-ft-long single-pass tube-in-shell heat ex- 
changer, which were expected to operate above 
1500°F, were constructed of niobium—1% zirconium 
alloy. The primary fluid, sodium or potassium, was 
electrically heated. Maximum oxygen partial pres- 
sure of ~10~° torr was maintained by baking out the 
system at 1400 to 2000°F with temperature cycling 
and by exposing the system to vacuum to ensure that 
all niobium components had completely outgassed. 

Boiling was difficult to achieve with this apparatus, 
and, when obtained, was unstable. Two distinct boil- 
ing regimes were observed: (1) a low-heat-flux mode 
with relatively large oscillations in flow, pressure 
drop, and fluid temperature with low (<10%) quality 
and considerable superheating; (2) a high-heat-flux 
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i 10° Btu/(hr)(sq ft)| mode that was associated with 
relatively high quality (> 35%) and little superheating. 
Although these two conditions were repeatable over a 
series of tests, no explanation of the reason for two 
coexistent flow regimes was evident. Apparently 
some source cavities for boiling that were dormant 
in the low-heat-flux regime were activated at ele- 
vated temperatures. 

The test series was conducted with shell Peclet 
numbers ranging from 50 to 350, and the overall 
heat-transfer coefficients for the boiler were found 
to be between those predicted for constant-tempera- 
ture and constant-heat-flux surfaces. The heat-trans- 
fer coefficient and vapor quality were obtained as 
functions of flow rate, and pressure-drop data were 
taken. Temperatures >1920°F were attained, and 
over 400 hr of running time at >1500°F were ac- 
cumulated. Flow rates were ~1000 lb/hr, and the 
temperature differential across the heat exchanger 
was ~ 250°F. Flow oscillations up to 50% were ob- 
served in the first boiling mode (probably amplified 
by the condenser). 


TURBOELECTRIC COMPONENTS 


Turboelectric design involves the development of 
lightweight machinery, usually turboalternators, op- 
erable at rotational speeds from 12,000 to 36,000 
rpm and at turbine inlet temperatures from 650 to 
1800°F. Temperature gradients of 600 to 700°F 
across a few inches are encountered, and close 
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Fig. 13 Typical two-phase boiling-potassium pressure-drop data.'” 
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clearances and manufacturing tolerances are man- 
datory. Material creep and coefficient-of-expansion 
data for liquid-metal turbine operation up to and 
beyond 10,000 hr at high temperatures are not abun- 


dant for the materials of interest; turbine design 
technique is pressing the limits of metallurgical 
technology. 


Turboallernatoy for SNAP-2. The model-V com- 
bined rotating unit (Table 2) developed by Thompson 
Ramo Wooldridge for the SNAP-2 system (Fig. 14) 
and the Mercury Rankine Power Conversion Program 
has a mercury-vapor-driven turbine, a pump, and an 
alternator on a common shaft.’ These rotating com- 
ponents of the mercury loop are also lubricated by 
mercury. Figure 15 shows the bearings for possible 
use in this system. 

The 4,.1-kw alternator is an 1800-cycle unit with a 
design efficiency of 87%. This three-phase device 
has a four-wire Y-connected circuit that is ex- 
ternally grounded neutral. 


The two-stage turbine consists of a_ subsonic 
impulse type first stage with 30% admission and a 
transonic impulse-reaction second stage with full 
admission. 


Two sets of impellers were tested for the cen- 
trifugal mercury pump. The first set has five vanes 
with a 10° inlet angle and a 31° discharge angle, 
and the second set has six vanes with a 45° dis- 
charge angle. Pump efficiency of 30% was obtained 
in tests at 45 lb/min flow. 
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Table 
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SPECIFICATIONS FOR THE COMBINED ROTATING UNIT 





Speed, rpm 

Turbine 
Inlet pressure, psia 
Inlet temperature, °F 
Discharge pressure, psia 
Discharge temperature, 


(min) 
Ib/min 


Efficiency, 
Flow rate, 
Alternator 

Frequency, cycles /sec 

Voltage, volts 
Line-to-neutral 
Line-to-line 

Output power, kw 

Power factor 

Duty cycle 


Efficiency at gross power, % 


IN SNAP-2 (REF. 5) 


1800 


63 + 9.5 
110 + 16.5 
4.1 

1.0 to 0.8 
Continuous 


87 


Mercury pump 
Inlet pressure to jets, 
psia (min) 5.9 
Inlet temperature, °F (max,) 
Discharge pressure, psia 525 + 25 
Flow, lb/min 50 +5 
Efficiency (with dynamic 
seal), % 41 
Bearings 
Power consumption, watts 
Alternator 
Turbine 
Thrust 100 
Flow rate, lb/min 
Alternator 8-10 
Turbine 8-10 
Thrust 





The combined rotating units have been subjected 
to design-verification tests, including endurance, 
shock, and vibration testing. Five units have ac- 
cumulated a total of > 30,000 hr of running time with 
more than 100 starts. Unit endurance-run times 
were 6637, 2558, 2049, 861, 3159, 1277, 760, 863, 
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Fig. 13 Mercury-lubricated bearings for the SNAP-2 com- 
bined rotating unit.* 


2500, 4767, and > 900 hr, for a total of ~ 25,000 hr. 
Degradation during the endurance runs has been very 
slight. 

The combined rotating unit has been run with 
preheating to 400°F with no detrimental effects noted; 
however, a preheating level of 200°F apparently 
resulted in some bearing wear. At near rated con- 
ditions, 4.2 kw(e) has been generated at a maximum 
mercury temperature of 1150°F. The unit has dem- 
onstrated its capability of operating for 2500 hr and 


Vol. 9, No. 4 


has passed all dynamic and performance tests, ex- 
cept that zero-gravity operation has not been verified. 


Turboalternator for SNAP-8. The SNAP-8 turbine 
is a four-stage axial-flow impulse machine with a 
rotational speed of 12,000 rpm. The first two stages 
have 38% admission, whereas the latter two stages 
have full admission. The turbine wheels are of 
Haynes Stellite 6B, and the inner housing is type 316 
stainless steel. Inlet temperature and pressure are 
1260°F and 265 psi; exit conditions are 265°F and 
16.5 psia. Turbine measured efficiency is 54% (the 
dry-gas efficiency is 58%). The turbine and alternator 
operate with separate but connected shafts. Overall 
efficiency is ~ 49.4% (Ref. 14). 

The turbine bearings are angle-contact preloaded 
units of iron—silicon—bronze, mounted in a back-to- 
back configuration. They are lubricated with 4P3E 
and include dynamic slingers so that flooded opera- 
tion is avoided. A portion of the turbine shaft is 
vented to space to prevent bearing lubricant from 
leaking into the mercury working fluid. 

The alternator is a one-piece straddle-mounted 
rotor with angle contact bearings lubricated with 
4P3E, which also cools the windings. This unit has 
been operated successfully up to a power level of 
53 kw(e) and has an efficiency of 87%. In testing, the 
turboalternator has run 830 hr and had 45 starts 
prior to the failure of a spiral-wound snap ring in 
one of the turbine-shaft labyrinth seals. The Stellite 
6B in the nozzle diaphragm and turbine wheel proved 
to be notch sensitive, and cracks were observed in 
these components. During operation the grain struc- 
ture of the Stellite 6B changed. Redesign eliminated 
all sharp corners and the separate blade shroud. A 
type 316 stainless-steel nozzle block was used, and 
Stellite 6B blades were inserted in the wheel. 

During test runs totaling 830 hr, the Stellite 6B 
blades did not erode, but performance for 10,000 hr 
has not been verified. Bearing-wear and space- 
sealing problems remain to be solved for the 10,000- 
hr duty cycle. Turbine critical speed is 20,000 rpm; 
however, the unit has been operated at 19,000 rpm 
with no significant increase in vibration. In general, 
performance requirements for this unit have been 
met, but endurance and reliability have yet to be 
achieved. 


Turbogeneralor for SNAP-50/SPUR. The turbogen- 
erator under development by AiResearch Mfg. Co. 
for the SNAP-50 system was intended to produce 
350 kw(e). The turbine design, Fig. 16, was for a six- 
stage 50% reaction device with a final stage diameter 
of 8.5 in. and a speed of 24,000 rpm. The final stage 
has an exit quality of 86% (Ref. 6). A single journal 
bearing is provided at each end of the turbine, and a 
thrust bearing is installed at the discharge end. 

A problem is the development of suitable bearings 
(Fig. 17). A single-stage Inconel X impulse turbine 
with supersonic nozzles and an efficiency of 50% was 
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Fig. 17 Dynamic bearing test rig uses a single-stage Inconel X impulse turbine.® 


employed in bearing tests. This wheel was lubricated 24,000 rpm. Table 3 summarizes the bearing-wear - 
by 550°F potassium, which was supplied through test experience. 

grooves in the bearing sleeves. Excellent bearing- 

test results were obtained with tungsten carbide A 200-Hp Two-Stage Potassium-Vapor Turbine. 
tilting-pad bearings in a type 316 stainless-steel Under a NASA contract General Electric Company 
housing. The tungsten bearings were run for over developed a 200-hp potassium turbine that repre- 
1000 hr continuously without evidence of wear, and, sented the third and fourth stages of a Rankine-cycle 
although subsequent restarts resulted in journal wear turbogenerator (Fig. 18) (Ref. 15). Turbine speed was 
(as a result of retaining-ring movement), the lubri- 19,200 rpm, and it operated with a pressure ratio 
cated bearings continued to operate smoothly at of 2.95. 


Table 3 FRICTION AND WEAR* OF BEARING MATERIALS IN 500°F POTASSIUM 


























Starting Running Celeulated 

Materials,# torque, torque, coefficient Weight 

rotor/shoe in.-lb av. in.-lb of friction change, mg Remarkst 
K94/Mo-0.5 Ti 18,3 13,1 0.266/0.190 +0.6/—4.0 Poor 
K601/Mo—0.5 Ti 17.4 10,9 0.253/0.159 -1.1/-3.7 Poor 
TZM/K94 > 26 Shutdown 0,380/ —40.0/—9.2 Very poor 
K94/TZM 17.0 10.9 0.247/0.159 +1.0/-—5.0 Poor 
K601/TZM 28.4 Fractured  0,410/ 
TZM/K94 26.2 14,0 0,381/0,.204 —59.0/—0.5 Poor 
Mo—0.5 Ti/K94 24.0 17.4 0.350/0.253 -—69.8/+0.2 Poor 
K601/TZM 28.4 Fractured 0,410 
TZM/W > 3 Shutdown 0.510 —108.4/—14.0 Very poor 
W/TZM 17.4 0.333/0.253 —15.4/-0.5 Poor 
K94/K94 15.3 13.0 0.222/0.189  -0.2/—0.4 Excellent 
K94/W 26.2 19.6 0.381/0.285  -666.3/—288.6 Poor 
W/TZM 1.3 4.0 0.062/0.058  —5.5/—0.7 Poor 
Mo-—0.5 Ti/W 2.1 0.03 —50.0/45.0 Poor 
K96/K96 8.3 2.75 0.121/0.025  -1.1/-0.1 Excellent 
K11/K11 4.36 2.18 0.063/0.052  —0.8/—0.3 E} 
K94/K94 4.36 : 0.063/0.063 +0.2/+0.2 E) 
K701/K701 2.18 ‘ 0.032/0.057 -—0.6/—0.2 E 
K701/K11 3.06 2.62 0.045/0.0389 -0.5/—0.2 E 

1 


K11/K701 4.36 3.92 0.063/0.057 —0.9/-—0. Excellent 











*50 in./sec; 0. <50 ppm. 

TTZM: Mo—0.5% Ti-0.1% Zr-—0.02% C; Mo-—0.5 Ti: Mo—0.5 Ti—0.02 C; W: pure 
are-cast W; Kll: WC—3 Co—0.1 TiC; K94: WC-12 Co—2(Ta, Cb)C; K96: WC-—12 
Co—2(Ta, Cb)C; K701; WC-—9.5 Co-0.1 TiC—4.0 Cr3Cy; and K601: 84.5 W—-10 Ta— 
5.5 C, no limit. 

tExcellent: Negligible wear, no material transfer, and polish both specimens, Poor: 
Excessive wear, surface roughened, material transfer, and high friction. Very poor: 
Test not completed due to excessive wear, friction, material roughening, and transfer. 
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Fig. 18 Saturated-vapor test turbine designed for 1600°F 


Of 


otassium at 92% inlet quality.” 
y 


One attempt was made to simulate the wet condi- 
tions of the latter turbine stages by spraying liquid 
into the inlet stream. This was not successful be- 
cause of the inability to obtain the fine droplets that 
would be present in an actual machine. 

Oil-lubricated bearings were employed during the 
tests to separate blade-erosion investigation from the 
equally important problem of bearing wear. Heat was 
supplied by a gas-fired boiler. Test conditions are 
shown in Table 4. 





Table 4 POTASSIUM—TURBINE TEST CONDITIONS'5 
Inlet temperature, °F 1450 —1580 
Turbine speed, rpm 15,400 —21,000 
Pressure ratio 1.4-—4.4 


Flow rate, lb/sec 1.6—2.1 
Efficiency, % 40-70 
Inlet quality, % 99 





Subsequent to the successful completion of per- 
formance tests, an endurance run was begun, and a 
total of 5000 hr (two runs) was logged. 

There was no visible second-stage blade erosion. 
A three-stage machine is being built for tests at 
higher last-stage moistures. 


Organic Working-Fluid System. A relatively low 
temperature Rankine turbine developed by Sundstrand 
uses Dowtherm as the working fluid.‘ This system 
consists of a boiler, turboalternator, regenerator, 
condenser, and pump in which the maximum fluid 
temperature attained is 653°F. The working-fluid 
serves as the lubricant for all rotating components. 

The shrouded, single-stage, Inconel X, radial- 
impulse turbine (Fig. 19) has a wheel radius of 3 in. 
and rotates at 24,000 rpm. It uses four supersonic 
(Mach 3) nozzles with a pressure ratio of 100 at an 
inlet pressure and temperature of 75 psi and ~ 650°F. 
Turbine efficiency of 54 to 61% has been attained. 
Because the expansion of the Dowtherm through the 
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turbine is a superheated process, no condensation 
problems are encountered. The turbine nozzles are 
set at a fixed tangential angle of 17°; power output 
is varied by changing nozzle size. 

A Pitot pump was employed since a minimum- 
leakage design was required for the very low (0.06 
Ib/sec) flows associated with this system. At 24,000 
rpm the pump exit pressure is 95 psi. Only partially 
successful, this design attained the specified pres- 
sure ratio but not the specified flow rate, a defi- 
ciency attributed to inlet cavitation effects. Net pump 
efficiency achieved was ~ 26%. 

The homopolar alternator consists of a main 
alternator and a parasitic unit that provides for load 
variation with constant turbine output. This 120- to 
280-volt three-phase 3200-cycle unit was initially 
intended to produce 1.5 kw(e) but later was re- 
designed to yield 4 kw(e). Design efficiency of the 
alternator is 85%. 

The bearings of the combined rotating unit, which 
consists of turbine, pump, and alternator, were 
tested with dry-start spin-ups to 24,000 rpm. Over 
137 starts and 905 hr of operation were accumulated 
on one set of bearings, and the feasibility of startup 
with dry bearings was verified. The turbine and alter- 
nator employed in these tests accumulated > 2200 hr 
of operation. 

The system was designed for flight use with a 
solar heat source. Therefore component weights are 








Fig. 19 Organic Rankine-cycleimpulse turbine wheel uses 
115 electrochemically machined Inconel X blades.‘ 
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minimal (Table 5). Further testing and modification 
of this system are in progress at Sundstrand. 

Follow-on programs have been initiated. One of 
the key problems in the long-term operation of an 
organic system is the stability of the working fluid 
with time and temperature. The AEC-sponsored 
programs at Sundstrand, Philco, and the National 
Bureau of Standards are investigating the long-term 
behavior of several candidate working fluids in simu- 
lated power-conversion loops. Under Air Force spon- 
sorship, Sundstrand is investigating the feasibility of 
a 2- to 10-kw(e) combined rotating unit. 


Table 5 


DOWTHERM RANKINE-CYCLE COMPONENT WEIGHTS 





Component Weight, lb 


Alternator housing and windings 32 
Alternator rotor 6 
Nozzle housing 7 
Pump housing 1.75 
Rear bearing housing 2 

1 


Pump and bearings 


Potal 51.25 


RADIATOR—CONDENSERS 


A number of designs have been tested and appear 
to be satisfactory. As with boilers, the limiting fac- 
tor is precise heat-transfer data. Radiator-tube 
vulnerability to meteors has also been studied briefly. 


Condensing Mercury in Low-Gravity Environ- 
ments. Two mercury-vapor-condensation test sys- 
tems have been set up at the NASA Lewis Research 
Center. The first of these (Fig. 20) was a nitrogen- 
cooled tapered-glass-tube apparatus that was tested 
in both 0- and 1-g environments. Glass tubes were 
selected to permit photography and to promote non- 
wetting flow for the tests. Test results included 4000 
frame/sec records of the condensing flow pattern. ' 










~ Cooli 


Condensing tube- 


Fig. 20 
Condensing mercury-vapor 
test rig, 


Pressure 
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manifold -/ 
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Pressure drops were correlated with the Martinelli 
parameter over flow rates from 0.028 to 0.049 Ib/ 
sec. The Martinelli correlation proved effective for 
high-quality flow (a low value of the two-phase flow 
modulus) but gave results significantly lower than the 
measured value for low-quality flow. The presence 
or absence of gravity had little effect on the mea- 
sured pressure drops. 

Tests were made with both wetting and nonwetting 
flow conditions (Fig. 21), and it was found that, at 
the lower flow rates, the pressure drop was only 
slightly affected by wetting. At higher flow rates the 
nonwetting pressure drop was lower than the wetted- 
tube pressure drop. 

A second series of tests was conducted with a 
NaK-cooled counterflow tapered-tube mercury sys- 
tem. As predicted, it was found that the mercury 
inlet static pressure initially dropped off rapidly 
with increasing condensing length but eventually 
leveled out. 


Mercury Rankine Radiator —Condenser. Atomics 
International conducted a series of mercury-radiator 
development tests for a 3- to 5-kw(e) Rankine system 
with 19.8 lb/min flow. Quality of the 625°F mercury 
entering the prototype radiator was 97%. An exit 
temperature of 371°F was required with a flow- 
passage pressure drop of 2.5 psi against a 9-psi 
backpressure.’ 

The test radiator is an array of tapered Haynes 25 
steel tubes with aluminum cooling fins. Versions with 
24 and 36 tubes were tested (preheat, operate, drain, 
and shutdown). The condensing length of the array 
was ~120 in. Although the radiator temperature dur- 
ing operational tests was nearly constant at 625°F 
over the initial 100 in. of tube length, a very rapid 
drop to the rated exit temperature occurred in the 
last 10 to 20 in. of tubing. The configuration tested 
appeared to be successful, but experiments with other 
configurations were planned. 
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Fig. 21 Flow configurations at interface locations for 
wetting mercury condensation in copper tube (upper) and 
nonwetting condensation in stainless-steel tube." 


Two full-scale horizontal radiator —condensers 
Simulating zero-gravity operation were tested. Over 
a range of conditions, no instabilities were noted. 
Satisfactory operation was also achieved against 
gravity by tilting the installation slightly. 


Mullitube Condenser for SNAP-8. Aerojet-General 
Nucleonics has tested a single-phase counterflow 
tube-in-shell NaK-—mercury device containing 73 
tapered 9% chromium—1% molybdenum steel tubes 
in a type 410 stainless-steel shell.'® The 40-in.-long 
tubes are 0.25 in. in diameter at the inlet end and 
0.20 in. at the exit. The total flow in the condenser 
is 11,500 lb/hr, and the variation in fluid quality is 
about linear with the tube length. The condenser 
pressure drop was successfully correlated with a 
modified Martinelli equation. 

The only significant fabrication problem encoun- 
tered in the assembly of the condenser was the weld- 
ing of the 4° tapered tubes to the tube header support 
plate. A tapered-backbraze technique proved to be 
successful. 

Tests showed that 420 kw(t) could be transferred, 
of which 15% was caused by subcooling of the con- 
densate. Multitube-array performance differed little 
from results obtained earlier with a single-tube 
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apparatus. This condenser appears to have satisfied 
all design requirements. 


Direcl-Condensing Potassium Radiators. A potas- 
sium, Single-fluid, radiator—condenser for a space 
power system was tested at ORNL.!® Each tube has 
copper radiator fins and is partially enclosed by a 
back reflector. A set of shutters attached to the re- 
flector can be employed to open or completely en- 
close each condenser tube. These shutters, closed 
for launching, can be kept shut for low-power opera- 
tion to control tube temperature and limit condensa- 
tion. For operation at rated power, the shutters are 
opened and form part of the reflector. 

Operational tests indicated satisfactory perfor- 
mance; no maldistribution or instability was ob- 
served. Data-gathering equipment included infrared 
films of the 800°F radiator surface. 


Polassium-Condensing Tests. A 50-kw(t) condens- 
ing heat-transfer loop (Fig. 22) has been operated 
at the General Eleciric facility in Evendale, Ohio.” 
Single tubes %, and * in. in inside diameter were 
flow tested with and without tapered-pin and helical 
inserts. 

Potassium temperatures ranged from 1100 to 
1400°F, and heat fluxes from 3 x 10‘ to 3 x 10° Btu/ 
(hr)(sq ft) were obtained. Flow velocities of 100 to 
1700 ft/sec were produced. 
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Fig. 22 Condensing potassium has been studied in this 
50-kw(t) loop? at heat fluxes up to 3 x 10° Btu/(hr)(sq fb. 








212 POWER REACTOR TECHNOLOGY 


The Martinelli correlation successfully predicted 
pressure variations, and no evidence of flow in- 
stability was encountered. The data taken (Fig. 23 
is typical) included the axial temperature distribu- 
tion, heat-transfer coefficients vs. temperature, the 
variation of the Martinelli parameter with tempera- 
ture, and the radial temperature variation in the 
condensing-tube walls. 
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Fig. 23 Potassium-condensing data obtained with the test 
loop shown in Fig. 22 (Ref. 20). 


A three-loop potassium-condensing system was 
assembled at NASA Lewis Research Center.”' This 
system consisted of a NaK primary loop, a two-phase 
potassium secondary loop, and a NaK cooling loop. 
Each NaK loop flowed 35,000 lb/hr, whereas the two- 
phase potassium loop design flow rate was 600 lb/hr. 

The potassium condenser was a type 316 stainless- 
steel tube-in-shell assembly 16 ft long. The tubes 
had an outside diameter of 0.5 in. anda wall thick- 
ness of 0.03 in. The inside diameter of the shell was 
2.0 in. Actual test conditions involved 1400°F potas- 
sium flow at 300 to 500 lb/hr with an inlet quality 
of 88%. Heat fluxes as high as 450,000 Btu/(hr) (sq ft) 
were observed. 

A 14- by 3-ft radiator —condenser was also tested 
during this study. The radiator, ~ 10 ft long, con- 
sisted of an array of nine parallel, finned tubes with 
a wall thickness of 0.042 in. and an inside diameter 
of 0.5 in. (Fig. 24). The radiator was mounted ina 
cooled vacuum chamber to provide an environmental 
pressure of 107 mm Hg. Heat fluxes ranging from 
23,000 to 49,000 Btu/(hr) (sq ft) were obtained. 
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Fig. 24 Potassium radiator—condenser with nine finned 
tubes?! 


A typical condensing radiator wall-temperature 
profile was obtained in these tests over most of the 
length of the assembly; subcooling occurred near 
the radiator exit. With a vapor inlet temperature 
of 1426°F and inlet quality of 83%, a flow of 569 lb/hr 
was obtained. Pressure drop, temperature, and heat 
flux were measured. Satisfactory operation of the 
radiator—condenser was obtained over a range of 
pressures that included test conditions with a near- 
zero pressure drop. 


Radiator Meteor Vulnerability. A program con- 
ducted by the NASA Lewis Research Center has 
attempted to provide some basic empirical data that 
would permit predicting how deeply a high-velocity 
particle would penetrate the wall of a radiator tube.” 

Pyrex projectiles were fired into flat-plate targets 
of aluminum, type 316 stainless steel, niobium, 
nickel, and cobalt and into tubes of 2024T6 aluminum 
and type 316 stainless steel. Penetrations were cor- 
related by an empirical relation. The flat-plate data 
were not, in most cases, very different from cor- 
responding tube data. However, no general pattern 
of penetration was evident in tests that covered a 
range of penetrations from mild dimpling through 
severe spalling and perforation. These modes of 
damage are illustrated in Fig. 25. 

In a series of impacts performed with lined tubes, 
spalling occurred with thick liners, but thin liners 
either dimpled or were penetrated. If mild dimpling 
could be tolerated, tube and liner weight could be 
reduced by about 45%. 
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A series of tests with beryllium armor over type 
315 stainless-steel tubes revealed that damage could 
be expected consistently with unreinforced beryllium. 
A reinforcing scheme wherein kinked wire or steel 
screening was embedded in the beryllium considera- 
bly increased the efficacy of the armor. For a 500- 
day mission in which dimpled but not spalled radiator 
tubes would be tolerated, the minimum radiator spe- 
cific weight [lb/kw(e)| would be obtained with lined 
beryllium tubes. In order of increasing weight, ATJ 
graphite, vanadium, and molybdenum follow beryl- 
lium. 


Damage mode: 


Closure Perforation Spall 








0.015 in. liner 0.065 in. liner 


0.028 in. liner A 





Inch 


Fig. 25 Limiting damage modes in type 316 stainless-steel 
tubes armored with cast aluminum. Tubes were at room 
temperature and were impacted by %/49-in.-diameter pyrex 
projectiles at 25,000 ft/sec (Ref. 22). 


Conclusions 


Most of the current work in the field is directed 
toward building a backlog of component experience. 
Because of the complexity and cost of this kind of 
power system, development must proceed along a 
rather narrow front and in a somewhat piecemeal 
fashion. The total effort reported in this article was 
under the aegis of only four government and seven 
commercial organizations. 

The prediction of heat-transfer coefficients in two- 
phase flow, wherein semiwetting conditions prevail, 
and the development of high-temperature turbines 
(1200 to 1800°F) and high-speed high-temperature 
bearings constitute the foremost technical problems. 
These represent a very straightforward projection 
of conventional rotating-machinery problems into a 
more severe and stringent environment. There is 
also the reliability problem—these Rankine sys- 
tems must operate for relatively long periods with- 
out maintenance. At present an engineering approach 


ENERGY CONVERSION 213 


to the solution of these problems is, for better or 
worse, being followed. Perhaps the seeming under- 
emphasis on basic heat-transfer and dynamics re- 
search is only an apparent one, and fundamental 
design principles as well as practical techniques 
may emerge from current Rankine-cycle work. 
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